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INTRODUCTION 


"FLYIT3"  is  a 3-degree-of-f reedom  trajectory  simulation  program 
written  in  FORTRAN  IV.  The  three  degrees  of  freedom  are  range,  altitude, 
and  pitch  attitude.  Pitch  attitude  (or  angle-of-attack)  is  a pseudo 
degree  of  freedom,  in  that  its  rate  of  change  is  unaffected  by  the  mass 
properties  of  the  vehicle.  In  this  respect,  the  vehicle  is  treated  as 
a particle,  the  forces  acting  on  which  are  controlled  by  its  attitude, 
among  other  things. 

The  program  is  useful  in  simulating  vertical-plane  trajectories 
involving  take-offs,  climbs,  cruises,  dives,  landings,  and  rocket  boosts. 
The  program  assumes  that  the  earth’s  surface  is  "flat,"  in  that  no  force 
is  applied  to  the  vehicle  (other  than  gravity)  to  prevent  growth  of  the 
altitude  coordinate  as  the  range  increases.  The  program  is  particularly 
well  suited  to  the  synthesis  and  analysis  of  air-to-surface  and  surface- 
to-surface  missiles  (especially  cruise  missiles)  and  aircraft  trajectories. 
FLYIT3  should  be  used  for  preliminary  study  only.  It  is  not  suitable 
(due  to  its  ignorance  of  rigid-body  effects,  primarily)  for  air-to-air 
trajectories,  terminal  engagements,  2-body  problems,  or  altitudes  over 
100,000  feet  MSL.  This  manual  is  written,  in  part,  to  familiarize  the 
user  with  FLYIT3's  limitations  as  well  as  its  capabilities. 

FLYIT3  is  easy  and  inexpensive  to  use,  is  compatible  with  batch  and 
demand  terminals,  and  at  NWC  can  automatically  produce  a plot  of  any 
trajectory.  The  program  is  simply  structured,  permitting  a user  to 
expand,  contract  or  alter  its  capabilities.  In  most  cases,  the  user 
must,  in  fact,  provide  his  own  coding  to  describe  his  vehicle's  aero- 
dynamic and  propulsion  characteristics.  However,  this  manual  provides 
comprehensive  guidance  in  this  area.  No  coding  need  be  provided  by  the 
user  to  "control"  the  trajectory.  Such  control  is  provided  by  careful 
selection  of  the  input  data.  All  known  "bugs"  in  the  program  have  been 
eliminated  (although  it  is  not  guaranted  that  there  are  none).  In  most 
cases  of  difficulty,  the  content  or  format  of  the  input  data  is  at  fault. 

This  manual  assumes  the  reader's  familiarity  with  FORTRAN  and  the 
UNIVAC  1110  executive  system.  Operation  is  described  in  terms  of  demand 
terminal  usage,  though  batch  mode  is  no  less  compatible. 
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ACCESS  TO  THE  PROGRAM 


At  NWC,  FLYIT3  is  stored  in  the  read-only  file  FLYIT*3.  The  program 
is  also  stored  in  a file  of  the  author  and  as  a card  deck.  Users  must 
copy  FLYIT*3  into  their  own  files  (permanent  or  temporary)  before  use. 

A minimum  of  200  tracks  should  be  assigned  to  the  file  to  permit  addi- 
tions or  alterations  to  be  made  to  the  program  in  the  user's  file.  In 
case  of  difficulty,  contact  the  author  or  inquire  of  the  Propulsion 
Analysis  Branch,  Code  3241. 

STRUCTURE  OF  FLY  IT 3 


The  source  program  is  written  in  FORTRAN  IV.  Execution  and  flow-of- 
control  is  performed  by  the  main  program  element,  FLYIT3.  The  main  pro- 
gram calls  subroutines  from  other  functional  groups  as  needed  to  supply 
data  with  respect  to  aerodynamics,  propulsion,  etc.  The  functional 
groups  are:  atmosphere,  trajectory  controls,  aerodynamics,  rocket 
boosters,  and  sustain  propulsion.  Each  of  these  groups  contains  several 
subroutines,  one  or  more  of  which  is  generally  used  for  each  trajectory. 
Each  subroutine  is  (has  been)  compiled  into  a binary  relocatable  element 
with  the  same  name  as  its  symbolic  source  element. 

The  file  FLYIT*3  contains  a runstream  element  titled  FLYIT*3. 
REASSEMBLE.  As  the  name  suggests,  this  element  is  @ADDed  to  MAP  the 
relocatable  elements  into  an  executable  element  named  EXECUTABLE. 
REASSEMBLE  need  not  be  used  unless  a source  program  element  has  been 
altered  and  recomplied  since  the  last  run. 

Another  runstream  element,  DATA,  is  @ADDed  to  execute  the  program. 
DATA  contains  the  data  describing  the  vehicle,  the  trajectory  to  be 
flown,  and  the  output  format.  In  addition,  DATA  contains  instructions 
to  describe  the  trajectory  plots  to  be  produced,  if  desired.  FLYIT3 
uses  NWC's  SHOMAT  plotting  program  to  produce  plots  at  minimum  cost. 

See  sections  titled  "DATA"  and  "PLOTTING." 

Each  functional  subroutine  group  is  described  in  some  detail  in  the 
following  sections. 


MAIN  PROGRAM:  FLY IT 3 


Subroutine  FLYIT3  is  listed  in  Appendix  A.  It  begins  with  condensed 
input  instructions  to  the  user  in  the  form  of  comment  lines.  Upon 
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execution  (0ADD  FLYIT*3. DATA)  this  subroutine  informs  the  user  it  is 
ready  to  read  data  by  printing  "enter  data."  The  data  is  then  read  (3 
lines)  from  the  runstream  element  DATA  without  further  interaction  with 
the  user.  FLYIT3  then  establishes  the  output  format  based  on  the  data 
received.  Initial  conditions  are  established,  also  based  in  part  on  the 
input  data.  When  certain  output  options  are  specified  in  the  data, 

FLYIT3  prints  the  input  and  output  data.  Input  options  are  discussed 
under  "DATA"  below.  Next,  FLYIT3  enters  four  nested  DO  loops  within 
which  the  trajectory  is  simulated.  The  outer  loop  in  effect  defines  a 
limit  on  the  number  of  "pages"  of  output  printed,  via  a "max  time"  limit. 
The  next  loop  limits  the  number  of  output  lines  printed  to  46,  at  which 
time  new  output  column  headings  are  printed  and  output  printing  may 
continue.  On  demand  terminals,  this  results  in  page-length  segments  of 
output  suitable  for  binding.  On  batch  printers,  the  new  output  headings 
are  printed  at  the  top  of  a new  page.  The  next  loop  controls  the  print 
interval  in  terms  of  the  number  of  iterations  calculated  before  a line 
of  output  is  printed.  The  innermost  loop  controls  the  frequency  with 
which  new  atmospheric  properties  are  calculated.  (This  feature  saves 
time  and  cost  under  conditions  where  atmospheric  properties  change 
slowly  or  not  at  all.) 

FLYIT3's  next  move  (within  the  innermost  loop)  is  to  call  one  of 
three  atmosphere  subroutines,  ATMSTD,  ATMHOT,  or  ATMCLD  (if  required). 
These  subroutines  supply  to  FLYIT3  the  properties  of  standard,  hot-day, 
and  cold-day  atmospheres. 

The  next  step  is  to  call  one  of  Che  subroutines  Ml,  M2,  M3...M10, 
which  selects  the  vehicle  angle-of-attack  required  to  fly  the  desired 
trajectory.  In  some  cases,  the  M-subroutines  in  addition  constrain 
accelerations,  flight  path  angles,  altitudes,  etc.  when  required.  The 
subroutine  called  in  this  group  depends  on  the  value  of  the  variable 
MGO.  MGO  is  always  equal  to  the  subscripted  variable  M(MM) . The  sub- 
script MM  is  the  position  in  a data  string  of  the  desired  sequence  of 
trajectory  segments.  Each  value  of  M()  from  1 thru  10  in  the  data  string 
corresponds  to  a particular  mode  of  flight  (climb,  dive,  cruise,  etc.). 
The  modes  indicated  by  M()  are  briefly  described  below: 

M*1  Runway  take-off  at  max  power  from  field  elevation  XLH. 

M»2  Air  launch  at  XLH  altitude,  XLMN  Mach,  XLFEE  flight  path  angle, 
and  climb  at  CLMN  Mach  DFEEC  flight  path  angle  and  continuously 
reduce  flight  path  angle  if  necessary  to  maintain  speed. 

M*3  Cruise  at  DCH  altitude  and  CMN  Mach  until  RC  range  is  reached 
or  fuel  is  exhausted. 

M“4  Dive  at  -DFEEDD  flight  path  angle  until  pull-up  is  initiated. 

M”5  Pull-out  from  dive  at  GPU  load  factor  and  cruise  at  DTH 
altitude  and  DMN  mach. 

M«6  Fly  at  constant  velocity  until  DCH  altitude  or  RC  range  is 
reached  or  impact. 

M“7  Launch  with  flight  path  angle  DFEEC  first  100  feet. 

M=3  Fly  at  body  angle  BAC  or  BAD  degrees  until  DCH  is  reached. 
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M=9  Approach  Co  landing  at  glide  slope  DFEEDD,  runway  elevation  DTH, 
approach  speed  DMN  Mach. 

M*10  Directs  the  vehicle  to  fly  along  the  radar  horizon  when  the 
horizon  is  intercepted  in  M-3  cruise.  Target  antenna  is 
assumed  to  be  at  100  feet  above  sea  level.  Range  to  target 
must  be  input  as  RC  (nmi) . 

If  used  at  all,  M=1  and  M=7  must  be  first  in  the  sequence. 

If  used  at  all,  M=9  will  be  last  in  the  sequence. 

No  M value  may  be  used  more  than  once  in  a trajectory  except  M*2, 

3,  4,  6,  8. 

M=4  must  be  used  in  the  sequence  45,  or  be  the  last  in  the  sequence. 

M*10  if  used  at  all  must  be  last  in  the  sequence. 

(The  above  was  excerpted  from  the  input  instructions  contained  in  FI.YIT3.) 
A typical  data  string  (included  on  the  first  of  the  three  lines  of  input 
data)  might  be  1834560000.  In  this  case,  M(l)=l,  M(2)=8,  M(3)=*3,  etc.  ; 
other  input  data  indicates  cruise  altitudes,  climb  and  dive  angles,  Mach 
numbers,  etc.  for  each  segment  (M  value)  in  the  trajectory.  Input  data 
is  discussed  in  detail  under  "DATA"  below. 

Upon  receiving  commands  from  the  trajectory  controls  group,  FLYIT3 
calls  for  vehicle  aerodynamic  characteristics  from  the  aerodynamics  group, 
subroutines  AEROl,  AER02...  AER015.  The  subroutine  called  is  indicated 
by  the  input  data  value  NAERO.  For  example,  if  NAERO  is  12,  AER012  will 
be  called.  The  user  is  responsible  for  the  contents  of  AER012,  or  which- 
ever subroutine  he  uses.  Fifteen  AERO  subroutines  are  provided  so  that 
each  user's  file  may  constitute  a library  of  aerodynamic  characteristics, 
any  of  which  may  be  selected  by  a simple  change  in  the  data  value  NAERO. 
The  same  discussion  applies  to  the  boost  propulsion  group  and  sustain 
propulsion  group,  discussed  below. 

FLYIT3's  next  call  is  to  the  booster  propulsion  group,  subroutines 
BSTR1  thru  BSTR15.  The  controlling  input  data  value  is  NBSTR.  If  NBSTR 
is  blank,  zero  or  negative,  a booster  is  assumed  to  be  absent  throughout 
the  trajectory.  If  a booster  is  used,  the  user  is  responsible  for  supply- 
ing its  thrust-time  profile,  specific  impulse,  burnt  weight  (which  will 
be  automatically  "dropped"  at  the  end-of -boost) , and  incremental  aero- 
dynamics, within  the  subroutine  he  uses. 

FLYIT3  then  calls  for  sustain  propulsion  information  from  subroutines 
SUST1  thru  SUST15.  Again,  the  user  is  responsible  for  the  contents  and 
selection  of  these  subroutines.  The  controlling  input  data  value  is 
NSUST.  Twenty  subroutines  are  available.  If  no  sustainer  is  present, 
NSUST  must  be  blank,  zero,  or  negative.  If  a sustainer  is  used,  the 
input  data  value  WO  (launch  weight)  must  be  greater  than  the  input  value 
WE  (empty  weight).  The  difference  in  these  values  must  include  sustain 
fuel  or  propellant  and  booster  weight.  The  sustainer  subroutines  must 
set  sustain  thrust  equal  to  zero  when  the  vehicle  weight  W becomes  less 
than  the  empty  weight. 
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FIG.  1.  Full-list  Output  Data  Format. 
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FLYIT3,  with  all  the  information  thus  obtained  from  subroutine  calls 
each  iteration,  calculates  the  vehicle  accelerations,  velocities,  and 
positions.  If  the  altitude  becomes  negative  or  exceeds  105,100  feet  MSL, 
or  the  Mach  number  exceeds  5.6,  the  program  will  terminate.  Range  and 
altitude  values  are  stored  for  subsequent  plotting  and  output  data  are 
printed  at  selected  iteration  intervals.  The  iteration  time  interval  is 
normally  0.25  second,  but  is  reduced  during  boost  (0.025)  and  during 
periods  of  maneuvering,  and  lengthened  (to  0.5  second)  during  cruise. 
Finally,  when  the  trajectory  has  been  claculated  to  completion,  FLYIT3 
calls  the  SH0MAT  (Ref.  1)  plotting  program  if  plots  have  been  requested 
by  the  input  data. 


FULL-LIST  OUTPUT:  0UTPT1 


0UTPT1  is  a subroutine  used  to  print  a detailed  summary  of  the  input 
data  (Figure  1).  This  format  is  not  normally  used  for  routine  demand 
terminal  work.  The  alternative  output  format  for  listing  the  input  for 
each  trajectory  is  defined  and  printed  by: 

7234560000  10114500 

1400.  900.  1.00  15.  .47  1.00  5.0  .01  .00  .0 

25.0  1.50  15.0-10.0  20.0  1.70  100.  10.0  .90  .1  .90 

STANDARD  ATMOSPHERE 

FIG.  2.  Demand-Mode  Output  Data  Format. 


ATMOSPHERES 


ATMSTD 

Subroutine  ATMSTD  generates  air  temperatures  and  densities  in 
accordance  with  the  U.S.  Standard  Atmosphere,  1962.  The  densities 
calculated  are  in  error  less  than  0.1  percent  below  36,152  feet  MSL  and 
less  than  1.0  percent  from  36,152  feet  to  over  120,000  feet.  (See  Figure 
3).  The  temperatures  calculated  are  in  error  by  less  than  0.05  degree  R 
up  to  105,000  feet  MSL.  Pressure  is  calculated  from  density  and  tempera- 
ture in  subroutine  FLYIT3.  The  resulting  pressure  is  in  error  by  less 
than  0.9  percent  up  to  108,000  feet  MSL.  Atmospheric  properties  are 
most  accurate  in  the  troposphere,  where  most  trajectories  are  flown. 


ATMHOT 


This  subroutine  supplies  air  data  in  accordance  with  the  MIL-STD-210A 
Tropical  Atmosphere.  No  error  analysis  has  been  done  on  this  model,  but 
errors  should  not  exceed  2 percent  or  0.2  degree  R up  to  100,000  feet. 
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ATMCLD 

This  subroutine  supplies  air  data  in  accordance  with  the  MIL-STD-210A 
Polar  Atmosphere.  Densities  and  pressures  are  accurate  to  1.0  precent 
below  40,000  feet,  but  accuracy  diminishes  at  higher  altitudes.  Tem- 
peratures are  accurate  within  0.2  degree  R. 

ATMSTD,  ATMCLD,  and  ATMHOT  are  listed  in  Appendix  C. 


TRAJECTORY  CONTROLS 


The  Ml,  M2,  . . . M10  subroutines  control  the  trajectory  shape  in 
accordance  with  the  user's  data  entries,  and  are  thus  the  real  heart  of 
the  program.  The  logic  and  flow  of  control  within  these  subroutines  is 
quite  complex  and  the  user  is  cautioned  against  casual  tampering.  The 
user  is  advised  to  become  familiar  with  this  section  and  the  section 
titled  "DATA"  before  attempting  to  use  the  program.  It  is  surprising, 
even  to  the  author,  what  a wide  variety  of  trajectories  is  capable  of 
being  simulated  by  FLYIT3  with  judicious  choice  of  input  data.  Numerous 
examples  are  displayed  in  Appendix  I.  The  M-subroutines  are  listed  in 
Appendix  D. 


Ml 


This  subroutine  simulates  a runway  take-off.  Propulsion  can  be 
booster,  sustainer,  or  both  in  series  or  parallel.  The  runway  is  assumed 
to  be  level  and  its  elevation  above  sea  level  is  the  input  data  value  XLH, 
in  thousands  of  feet.  Liftoff  occurs  when  the  airspeed  reaches  1.20  times 
the  airspeed  at  which  lift  would  equal  weight  at  the  angle-of-attack  limit, 
specified  by  the  input  data  value  ALPL  (degrees).  At  the  time  of  liftoff, 
control  of  the  trajectory  is  transferred  to  the  next  M-subroutine  as 
indicated  in  the  input  data  sequence  M(l),  M(2) , . . .M(10) . 


M2 


The  M2  subroutine  is  used  for  climbing  and  diving  at  a specified 
flight  path  angle.  M2  can  be  used  as  the  first  in  the  M-sequence  to 
simulate  an  air- launch.  (M4,  M6,  M7,  M8,  and  M9  can  also  be  used  for 
air-launch.  For  any  air-launch,  the  launch  altitude  is  specified  by  the 
input  data  value  XLH,  the  launch  Mach  number  by  XLMN,  and  the  launch 
flight  path  angle  by  XLFEE.  "Air-Launch"  could  include  cases  in  which 
prior  to  "time  zero"  there  occurred  a boost  phase,  catapult  launch, 
ordnance  release,  etc.  In  this  way,  a complex  trajectory  can  be  pieced 
together.  The  capability  of  FLYIT3  rarely  requires  piecing,  however.) 
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The  flight  path  angle  at  which  the  climb  or  dive  is  to  proceed  is 
entered  as  the  input  data  value  DFEEC  (degrees  from  horizontal).  In 
addition,  a climb  Mach  number  can  be  entered  as  data  (CLMN) , which  sustain 
propulsion  can  be  throttled  to  maintain  (if  the  user's  propulsion  sub- 
routine is  so  programmed).  In  M2,  if  the  propulsion  is  unable  to  main- 
tain the  specified  climb  angle,  the  climb  angle  will  be  progressively 
lowered  so  as  to  maintain  approximately  constant  Mach  number  (to  prevent 
stall  or  "mushing"). 

M2  (also  M8)  anticipates  arrival  at  some  cruise  altitude  (DCH)  during 
the  climb.  The  vehicle  is  caused  to  initiate  a ballistic  (zero-lift)  turn- 
down at  approximately  the  time  which  permits  arrival  at  the  cruise  alti- 
tude at  a flight  path  angle  of  zero.  M2  transfers  control  to  the  next 
M-subroutine  in  the  M-sequence  when  fuel  is  exhausted,  range  exceeds  the 
input  data  value  RC  (nautical  miles),  or  cruise  altitude  is  reached.  If 
the  flight  path  angle  is  reduced  to  essentially  zero,  the  transition 
(usually  to  M3)  will  be  made  and  DCH  redefined  as  the  present  altitude. 


M3 


M3  provides  for  cruising  at  altitude  DCH  and  Mach  number  CMN . Stable 
altitude-hold  is  achieved  by  the  use  of  an  equation  which  represents  a 
proportional  feedback  control  system.  Two  error  signals  are  used,  alti- 
tude itself  and  the  flight  path  angle.  Damping  is  provided  by  mixing  the 
required  angle-of-attack  with  the  one  used  in  the  previous  iteration. 

The  user  is  cautioned  against  changing  the  gains  and  coefficients  in  the 
altitude-hold  equation,  as  they  have  been  optimized  with  respect  to  a 
broad  spectrum  of  vehicles  and  trajectories  over  many  years  of  use. 

The  M3  cruise  leg  of  a trajectory  is  normally  terminated  by  either 
of  two  events.  The  first  is  the  attainment  of  a range  specified  as  input 
data  value  RC  (nautical  miles).  The  second  is  fuel  exhaustion  (current 
vehicle  weight  less  than  vehicle  empty  weight  (WE)  pounds  in  data)  . When 
M3  is  terminated,  trajectory  control  is  transferred  to  the  next  M-subroutine 
in  the  M-sequence. 

M3  can  be  used  as  an  acceleration  leg  of  a trajectory,  if  that  leg 
is  to  be  followed  by  another  cruise  leg  at  a higher  altitude.  In  this 
case,  M3  terminates  when  the  Mach  number  CMN  is  attained  in  level  flight. 

The  second  M3  leg  will  be  flown  at  a Mach  number  of  CLMN  +0.5.  This 
feature  is  useful  when  a speed  much  higher  than  the  launch  speed  is 
required  for  good  climb  performance  to  a higher  altitude. 

If  M3  is  to  be  followed  by  an  M6  climb,  the  climb  will  be  initiated 
when  the  Mach  number  exceeds  CLMN.  M6  will  then  maintain  this  approxi- 
mate Mach  number  - see  under  M6. 
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M4  is  used  for  diving  at  constant  flight  path  angle.  M4  is  useful 
for  flying  a transition  from  a cruise  altitude  to  a lower  cruise  altitude 
(see  M5) , for  dives  to  impact  (...40...  in  the  M-sequence),  and  for  line- 
of-sight  trajectories  (e.g.,  WALLEYE).  The  diving  flight  path  angle  is 
entered  in  the  data  as  DFEEDD  (degees,  positive  for  dives).  Climbs  can 
also  be  performed  by  M4  by  entering  a negative  value  for  DFEEDD;  however, 
a smooth  transition  to  the  next  M-subroutine  in  the  sequence  cannot  be 
assumed.  (An  M4  climb  would  be  useful  for  simulating  a "zoom;"  for 
example,  just  prior  to  a parachute  deployment,  to  lower  the  airspeed.) 

The  input  data  value  DVMN  provides  for  a Mach  number  to  be  maintained  in 
the  dive  by  means  of  propulsion  throttling. 

Flight  path  angle  is  sought  and  maintained  (following  a ballistic 
turn-down)  by  a proportional  feedback  control  system  equation  which  uses 
flight  path  angle  as  the  error  signal  and  angle-of-attack  mixing  as  a 
damping  term.  In  dives,  M4  is  terminated  by  impact  (attainment  of  zero 
altitude)  or  by  arrival  at  a altitude  at  which  a pull-uD  must  be  initiated 
for  transition  to  an  M5  cruise  leg. 


M5 


M5  provides  for  a low  level  cruise  (dash)  leg  and  the  pull-up  thereto. 
Transition  to  M5  can  be  made  only  from  M4  (M5  is  also  used  by  M9  and  M10). 
The  altitude  for  this  dash  leg  is  input  as  DTH  (thousands  of  feet  MSL) , 
and  the  Mach  number  as  DMN.  If  DTH  has  the  same  value  as  DCH,  the 
M-sequence  ...345...  can  be  used  to  properly  transition  (smoothly)  from 
M3  to  M5  without  significant  disturbance  of  altitude.  This  sequence 
is  useful  for  flying  a constant-altitude  trajectory  with  a change  in 
speed  at  RC  range  from  launch. 

M5  provides  stable  altitude  hold  in  the  same  way  as  does  M3. 

The  pull-up  from  dive  to  level  flight  is  controlled  such  that  the 
maximum  normal  load  factor  is  GPU  g's,  the  angle-of-attack  limit  ALPL  is 
not  exceeded,  and  the  maneuver  is  initiated  (usually)  in  time  to  prevent 
substantial  undershoot  of  DTH.  Some  combinations  of  dive  angle,  speed, 
throttling,  aerodynamics,  and  DCH  cause  enough  undershoot  that  sea  level 
impact  threatens  the  trajectory.  In  this  case,  the  angle-of-attack  is 
increased  so  as  to  produce  load  factors  greater  than  GPU  in  the  interest 
of  salvaging  the  simulation.  To  avoid  this  problem,  it  is  recommended 
that  DCH  values  of  less  than  0.1  (100  feet  MSL)  not  be  used.  The  impact 
of  this  restriction  on  simulation  results  is  negligible,  considering  the 
typical  accuracies  represented  by  propulsion  and  aerodynamics  models. 
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M5  is  terminated  only  by  fuel  exhaustion,  at  which  time  a transition 
is  made  to  the  next  M-subroutine  in  the  M-sequence. 


M6 

The  M6  subroutine  flies  an  approximately  constant  speed  trajectory, 
holding  the  speed  at  which  transition  into  M6  is  made.  However,  the  speed 
is  maintained  by  continuous  adjustment  of  the  flight  path  angle  rather 
than  throttling  the  propulsion.  This  mode  of  flight  is  suitable  for 
launch,  climbs  and  dives,  and  cruise.  When  the  altitude  range  is  small, 
this  mode  simulates  constant  pitch  trim.  Transition  out  of  >16  can  be 
caused  by  a number  of  occurences,  such  as  range  exceeding  RC,  DCH  being 
reached,  etc. 


M7 


M7  provides  for  air-  or  surface- launch  in  which  the  direction  of 
motion  is  fixed  for  the  first  100  feet  of  travel  (rail  launch).  This 
option,  which  must  be  first  in  the  M-sequence  if  used,  is  convenient  for 
low-speed  or  zero-speed  launches  with  or  without  booster.  It  eliminates 
the  otherwise  troublesome  problems  of  launcher  tip-off  and  lack,  of  aero- 
dynamic control  and  lift  at  low  speeds. 

M8 


Subroutine  M8  provides  a means  of  climbing  and  diving  which  is  well 
suited  to  t; aversing  wide  spans  in  altitudes.  The  flight  mode  is  one  of 
a constant  pitch  attitude  with  respect  to  the  horizonal.  Two  pitch  angles 
are  input,  BAC  (body  angle  for  climbs)  and  BAD  (body  angle  for  dives). 

This  mode  of  flight  is  one  typically  used  by  a pilot  flying  an  aircraft, 
even  if  he  is  not  consciously  aware  of  using  it.  It  is  equivalent  to 
maintaining  a fixed  line  of  sight  between  the  pilot's  eyes,  some  ref- 
erence point  on  the  aircraft  (speck  on  the  windshield,  nose  highlight, 
etc.),  and  the  horizon.  For  fixed  thrust,  the  equivalent  airspeed  (EAS) 
will  be  approximately  constant.  Flight  at  constant  body  attitude  is 
relatively  easy  to  mechanize  in  an  airframe  control  system,  since  the 
only  sensor  required  is  a single-axis  gyroscope.  In  constant  attitude 
flight,  the  angle-of-attack,  airspeed,  lift,  drag,  and  flight  path  angle 
vary  in  an  inherently  stable  manner.  The  result,  in  an  atmosphere  whose 
density  decreases  with  altitude,  is  a tendency  for  the  flight  path  angle 
to  smoothly  approach  zero  and  the  speed  to  increase  as  high  altitudes 
are  approached.  The  tendency  for  flight  path  angle  to  increase  or  decrease 
depends  also  on  the  throttle  setting  and  the  actual  body  angle  selected. 

The  user  will  want  to  try  different  values  of  BAC  and  BAD  for  each  vehicle 
configuration  in  order  to  achieve  suitable  climb/dive  performance  (i.e., 
proper  variation  of  speed  with  altitude) . 
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M8  will  decide,  based  on  Che  alcitude  with  respect  to  DCH  at  the 
time  of  transition  to  M8,  whether  BAC  (climb)  or  BAD  (dive)  is  appro- 
priate. If  fuel  exhaustion  occurs  while  M8  is  active,  M8  will  fly  at 
BAD  to  impact.  Otherwise,  M8  will  fly  at  BAC  or  BAD  so  as  to  approach 
DCH.  However,  if  M3  has  previously  been  active  in  the  trajectory  or  an 
M2  climb  has  immediately  preceeded  current  M8  activity,  M8  will  seek  DTH. 
If  M3  is  the  next  indicated  subroutine  in  the  sequence,  but  the  range 
exceeds  RC  while  M8  is  active,  control  is  transferred  to  M3  at  that  time. 
(Subsequently,  M3  will  immediately  transfer  control  to  the  next  M-sub- 
routine  in  the  sequence.)  When  climbing  toward  DCH,  M8  will  anticipate 
arrival  and  initiate  a ballistic  turn-down  as  in  M2. 

M9 


M9  provides  for  simulating  a final  approach  to  a runway  landing  and 
the  subsequent  roll-out  with  braking.  The  runway  elevation  is  input  as 
DTH,  the  glide  path  angle  is  DFEEDD  (positive),  and  approach  Mach  number 
DVMN.  The  user  should  provide  suitable  aerodynamics  and  propulsion  for 
the  landing  configuration  (i.e,  flaps  and  wheels  down).  To  account  for 
braking  after  touchdown,  the  user  must  insert  a coefficient  of  rolling 
friction  (retarding  force/weight)  into  subroutine  FLYIT3.  The  variable 
name  for  this  coefficient  is  FRIC,  and  is  defined  as  a constant  (0.5) 
among  the  initial  conditions  of  FLYIT  (line  120,  approximately).  The 
user  must  edit  and  recompile  FLYIT3  whenever  FRIC  is  redefined. 

M9  will  print  a line  of  output  indicating  the  time  of  touchdown  and 
the  touchdown  flight  path  angle.  This  angle  is  different  from  the  glide 
path  angle,  because  a 1.1-g  flare  occurs  prior  to  touchdown.  Another 
output  line  is  printed  when  the  vehicle  rolls  to  a stop,  indicating  the 
total  runway  length  used. 

When  M9  is  indicated  in  the  M-sequence,  it  must  be  followed  by  a 
zero,  and  must  not  be  preceeded  by  45. 

M10 


Subroutine  M10  causes  the  vehicle  to  follow  a path  coincident  with 
the  radar  horizon  with  respect  to  a target.  The  range  from  launch  point 
to  the  target  must  be  input  as  RC  (nautical  miles).  The  target’s  radar 
antenna  is  assumed  to  be  100  feet  above  sea  level.  The  transition  to 
M10  must  occur  from  M3  cruise. 

The  radar  horizon  path  is  tangent  to  the  sea  surface  12.38  NM  from 
the  target.  Since  the  vehicle  cannot  be  expected  to  fly  at  -ero  altitude, 
it  is  caused  to  level  off  and  cruise  at  100  feet  MSL  when  the  altitude 
drops  to  100  feet.  At  this  time  M10  redefines  DTH  = 100  (feet),  and 
causes  a transfer  of  control  to  M5.  Depending  on  the  user's  propulsion 
coding,  the  vehicle  can  be  flown  at  Mach  numbers  of  DVMN  or  CMN  in  the 
descent,  and  DVMN,  CMN,  or  DMN  (among  others)  in  the  final  cruise  leg. 
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The  trajectory  will  terminate  in  M10  when  the  range  to  target  becomes 
zero.  M10  will  print  the  time-to-target  and  range  flown  at  impact. 


AERODYNAMICS:  AEROl  THRU  AERO 15 


Subroutines  AEROl  thru  AER015  are  available  to  the  user  to  define 
his  vehicles'  aerodynamic  characteristics.  The  file  FLYIT*3  includes 
AEROl,  AER02,  and  AER03  as  real,  text-contained  elements  which  serve  as 
examples  of  satisfactory  coding,  in  addition  to  possibly  being  useful 
to  the  user.  (See  Appendix  E for  listings.) 

Each  AERO  subroutine  must  provide  values  for  the  following  variables 
in  the  common  block: 

CDO  - zero-lift  drag  coefficient, 

CDI  - induced  (or  drag-due-to-lif t)  drag  coefficient, 

CL  - lift  coefficient  (not  normal-force  coefficient),  and 

CLAR  - the  rate  of  change  of  CL  with  angle  of  attack  in  radians 
(lift-curve  slope). 

All  these  coefficients  are  to  be  based  on  the  input  data  value  AREF,  the 
aerodynamic  reference  area  in  square  feet. 

In  order  to  calculate  these  variables,  the  AERO  subroutines  may  use 
any  of  the  variables  in  the  common  block  of  FLYIT3.  Most  typically,  these 
are: 


M(MM)  - the  active  M-subroutine  (or  MGO) , 

W - vehicle  weight, 

RO  - atmospheric  density, 

ALPHA  - angle-of-attack  in  degrees, 

ALPL  - angle-of-attack  limit  in  degrees, 

V - true  airspeed  in  feet  per  second, 

Q - dynamic  pressure  in  psf, 

H - altitude  in  feet  MSL, 

AMACH  - Mach  number,  and 

AC  - inlet  capture  area  in  square  feet. 

The  basic  aerodynamics  can  be  modified  by  the  user's  propulsion  sub- 
routines. This  is  usually  the  most  convenient  way  to  increment  zero-lift 
drag  to  account  for  detachable  boosters,  inlet  spillage  drag,  and  engine 
cold-flow  drag. 


* 
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Subroutine  FLYIT3  adds  the  basic  AERO-produced  and  incremental 
characteristics  and  calculates  the  lift  and  drag  forces. 


ROCKET  BOOSTERS:  BSTR1  THRU  BSTR15 


Subroutines  BSTR1  thru  BSTR15  are  available  to  define  the  character- 
istics of  rocket  boosters  which  may  be  used  on  a vehicle.  The  file 
FLYIT*3  includes  BSTR1,  BSTR2,  and  BSTR3  as  real,  text-containing, 
unclassified,  elements  which  serve  as  examples  of  satisfactory  coding, 
in  addition  to  possibly  being  useful  to  the  user.  (See  Appendix  F for 
listings.)  Many  booster  motors  have  been  modelled  by  and  are  available 
from  Code  3241. 

Each  BSTR  subroutine  must  provide  values  for  the  following  variables 
in  the  common  block: 

THRB  - booster  gross  thrust  in  pounds  (assumed  to  be  aligned 

with  the  angle-of-attack  datum,  usually  the  fuselage  axis), 

FISPB  - propellant  specific  impulse  in  lbf-sec/lbm, 

DT  - iteration  time  interval  during  boost  (this  value  over- 

rides values  established  by  the  M-subroutines  - 0.025  is 
recommended) , 

EWB  - the  weight  dropped  at  end-of-boost  in  pounds, 

Any  other  variable  in  the  common  block  can  be  modified  by  or  used 
by  BSTR.  These  typically  include: 

TIME  - time  accumulated  since  launch  in  seconds, 

CDO  - zero-lift  drag  coefficient  (usually  incremented  by  the 

additional  drag  of  the  booster) , 

H - altitude  in  feet  MSL  (used  to  augment  booster  thrust  and 

Isp  with  altitude) . 

The  input  data  value  NB  permits  multiple  boosters  to  be  used  on  one 
vehicle,  on  the  assumption  that  all  boosters  are  ignited  simultaneously. 


SUSTAINER  PROPULSION:  SUST1  THRU  SUST20 


Subroutines  SUST1  thru  SUST20  are  available  to  define  the 
characteristics  of  sustainers.  The  file  FLYIT*3  includes  SUST1  thru 
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SUST4  as  real,  text-containing , unclassified  elements  which  serve  as 
examples  of  satisfactory  coding,  in  addition  to  possibly  being  useful  to 
the  user.  (See  Appendix  G for  listings.)  Other  sustainer  engines  have 
been  modelled  and  are  available  from  Code  3241. 

Each  SUST  subroutine  must  provide  values  for  the  following  variables 
in  the  common  block: 

THRS  - sustainer  thrust  in  pounds 

FISPS  - sustainer  propellant  specific  impulse 

CDAOFF  - additional  zero-lift  drag  coefficient  resulting  from 

an  engine-off  condition  (must  be  referenced  to  AREF) 

Any  other  variable  in  the  common  block  can  be  modified  by  or  used 
by  SUST.  These  typically  include: 

H - altitude, 

AMACH  - Mach  number, 

3 

RO  - atmospheric  density  in  slugs/ft  , 

ALPHA  - angle  of  attack, 

Q - dynamic  pressure  in  psf, 

DT  - time  increment, 

TM  - sustainer  thrust  multiplier  (see  below) , 

2 

AC  - inlet  capture  area  for  unit  value  of  TM  in  ft“  (see  below), 

CMN  - desired  cruise  Mach  number, 

CLMN  - desired  climb  Mach  number, 

DVMN  - desired  dive  Mach  number, 

DMN  - desired  dash  Mach  number, 

PT  - freestream  total  pressure  in  psia, 

TT  - freestream  total  temperature  in  °R, 

AMACH1  - flight  Mach  number  prior  to  last  iteration, 

THRB  - booster  thrust, 

T - ambient  air  temperature  in  °R,  and 

MGO  - M-subroutine  currently  active. 

The  input  variable  TM  is  a sustainer  thrust  multiplier.  This  vari- 
able permits  the  sustainer  and  its  inlet  to  be  conveniently  scaled  in 
thrust.  If  TM  = 1.0,  the  engine  will  produce  the  performance  coded  in 
SUST.  If  TM  is  different  from  unity,  the  thrust  and  inlet  capture  area 
will  be  scaled  by  TM. 
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Each  sustainer  subroutine  must  contain  the  performance  of  an  inlet, 
if  applicable.  Inlet  additive  and  spillage  drag  is  usually  added  to  CDO 
in  SUST;  however,  the  variable  CDAOFF  is  available  as  an  incremental 
zero-lift  drag  coefficient  which  is  added  to  CDO  in  subroutine  FLYIT3. 

The  variables  ASPEED  and  TIT  (for  gas  turbines,  rotor  speed  and 
turbine  inlet  temperature,  respectively)  may  be  used  as  diagnostic  or 
informative  parameters.  They  are  printed  as  output  data  in  the  batch 
(full  list)  output  format. 

The  SUST  subroutines  are  responsible  for  setting  THRS  = 0.0  if  and 
when  the  vehicle  weight  W becomes  less  than  the  empty  weight  WE. 

Rocket  motors  may  be  modelled  as  sustainers  as  well  as  boosters; 
however,  the  coding  is  not  directly  interchangeable  between  BSTR  and  SUST. 

Contractor-supplied  engine  performance  coding  can  be  integrated  into 
or  called  from  a SUST  subroutine,  but  proper  interface  coding  must  be 
provided. 

SUST  is  responsible  for  engine  throttling  if  this  feature  is  desired. 
The  input  variables  CMN,  CLMN,  DMN,  and  DVMN  in  conjunction  with  MGO, 
AMACH,  and  AMACH1  indicate  whether  an  increase  or  decrease  in  thrust  is 
desired.  In  simple  throttling  equations,  it  is  recommended  that  a pro- 
portional control  system  with  damping  be  used.  The  quantity  (DMN-AMACH) , 
for  example,  would  be  the  error  signal  and  the  rate-of-change  of  Mach 
number  ( (AMACH-AMACH1) /DT)  would  stabilize  the  system.  To  approximate 
the  response  time  of  a real  sustainer  engine,  a limit  could  be  imposed 
on  the  rate  of  change  of  a critical  engine  parameter  (e.g.,  rotor  speed). 
The  throttling  equations  can  range  from  no  throttling  at  all  to  a full 
digital  simulation  of  a complex  fuel  control  system. 

The  author  disclaims  responsibility  for  the  accuracy  of  models 
contained  in  FLYIT*3  representing  real  engines,  existing  or  proposed. 

Code  3241  has  available  coded  models  of  many  engines. 


REASSEMBLE 


Element  REASSEMBLE  is  a runstream  element  which  when  @ADDed  maps 
the  relocatable  elements  in  FLYIT*3  and  the  library  file  PLT*PLOT&LIBRARY 
into  an  executable  program  element  named  EXECUTABLE.  REASSEMBLE  should 
be  @ADDed  prior  to  the  user's  first  run  attempt  and  after  any  source 
elements  in  FLYIT*3  have  been  edited  and  recompiled  (except  for  element 
DATA,  which  is  a runstream  element  and  is  not  complied).  FLYIT*3  should 
be  @PACKed  following  each  new  MAPping,  to  prevent  a storage  overflow. 
REASSEMBLE  is  listed  in  Appendix  H. 


DATA 


Element  DATA  is  a runstream  element  which  causes  the  absolute  pro- 
gram EXECUTABLE  to  be  executed,  and  which  contains  the  input  data  for 
the  run.  DATA  also  contains  plotting  instructions  for  the  SHOMAT  plotting 
program.  A successful  example  of  DATA  is  listed  in  Figure  4.  Line 
numbers  (1  thru  18)  have  been  included  for  reference,  but  are  not  part 
of  the  text  of  DATA.  Line  1 causes  EXECUTABLE  to  be  excuted.  Lines  2, 

3,  and  4 constitute  the  input  data  itself  (see  discussion  below).  Lines 
5 thru  17  constitute  instructions  for  SHOMAT  plotting.  If  no  plotting  is 
desired,  lines  following  line  5 are  ignored  and  can  be  omitted.  Line  18 
(i3E0F)  is  necessary  only  when  plotting  instructions  are  used.  It  signals 
the  end  of  the  plotting  instructions  and  causes  the  plotting  to  be  executed 
one  trajectory  per  page.  (9E0F  replaces  3FIN  as  discussed  in  Reference  1. ) 
The  user  is  referred  to  Reference  1 for  detailed  instructions  on  the  use 
of  SHOMAT. 


li@XQT  FLY  I T*3. EXECUTABLE 

2:723456  l_  _ JQ I 1 I 5_  0 0 

3:1400.1100.  1.00  15.0  .47  1.00  5.0  .01  0.9  45. 

4i  10.  1.5  15.  -10.  20.  1.7  150.  10.  0.9  0.1  ^0.9 

5* I DINFO  3241  G 3URGNER  PHONE  7301 
6:1 GSEWB 

7:T  L TRAJECTORY  SHAPE 

8»T  USERS_MANUAL_  EXAMPLES. 

9:  A RANGE,  NM 

1 0 :()  ALTITUDE,  FT  __ 

1 1 :0BJECT  2. , 1 .5, 9. 5, 6. 5 
12*SCALEX  0,300,  10,5,5  _ __ 

1 3 * T I CKX 

14ISCALEY  0,50000, 5000*2  ,.2 _ 

I 5 :TICKY 
I 6 : J0I N 
I 7 : END 

I3*GE0F  _ 


FIGURE  4 


A Successful  Version  of  Element  DATA. 
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INPUT  DATA  (LINES  2,  3,  AND  4) 

Lines  2,  3,  and  4 are  reprinted  and  dissected  in  Figure  5 for 
illustration.  Refer  to  the  section  on  trajectory  controls  for  supple- 
mentary descriptions  of  the  functions  of  each  variable.  When  constructing 
or  editing  lines  of  data,  pay  strict  attention  to  the  field  allotments 
(spacing) . Floating  point  data  should  have  decimals  punched.  Line  2 
should  concain  no  decimal  points.  In  the  short-form  output  data  format 
such  as  appears  in  Figure  2,  and  in  the  examples  of  Appendix  I,  DATA 
lines  1,  2,  and  3 are  repeated  as  output  in  a similar  format. 

Line  2:  The  format  is  (1011,815) . The  first  ten  fields  are  the 
M-sequence  which  indicates  qualitatively  the  trajectory  shape  and  modes 
of  flight.  (Lines  3 and  4 contain  values  which  quantify  the  trajectory 
and  the  vehicle.)  The  M-sequence  is  discussed  in  detail  under  "Trajectory 
Controls"  above.  The  next  field  contains  the  variable  NB,  the  number  of 
booster  motors  used  in  parallel.  This  number  is  meaningless  unless  a 
booster  unit  is  specified  by  the  variable  NBSTR.  The  next  field  contains 
the  variable  NATM.  A blank  or  zero  specifies  the  use  of  the  ATMSTD 
atmosphere  model  for  the  run,  "1"  calls  for  Topical  Atmosphere  (ATMHOT) , 
and  "2"  calls  for  Polar  Atmosphere  (ATMCLD) . The  next  field  contains  the 
variable  NAERO,  which  is  the  aerodynamics  subroutine  to  be  used;  i.e. , 

"1"  calls  AER01,  "2"  calls  AER02,  etc.  The  next  field  contains  NBSTR. 

"1"  calls  subroutine  BSTR1,  "2"  calls  BSTR2,  etc.  The  next  field 
contains  NSUST.  "1"  calls  SUST1,  "2"  calls  SUST2 , etc.  The  next  field 
contains  the  variable  NOUTPT.  The  allowable  options  are  0,  1,  2,  3,  4, 
and  5.  "0"  or  blank  specifies  the  132-column  output  format  (Figure  6) 

(batch-mode)  with  a short-form  listing  (Figure  2)  of  the  input  data. 

"1"  specifies  batch-mode  output  format  (Figure  5)  with  input  data  fully 
described  by  the  full-list  format  (Figure  1;.  "2"  specifies  option  "1" 

plus  a plot  of  the  trajectory.  "3"  specifies  the  70-column  demand-mode 
output  format  with  a full-list  description  of  input  data.  "4"  specifies 
the  short-form  demand-mode  output  format  (Appendix  I) . "5"  specifies 

option  "4"  plus  a plot.  The  next  field  contains  NPRT.  If  a positive 
number  is  placed  in  the  field,  the  trajectory  output  data  will  be  printed 
four  times  mofe  frequently  than  usual,  allowing  a "closer  look"  at 
maneuvers,  boost  phase,  etc.  The  last  field  contains  NFEET.  If  a 
positive  number  is  in  this  field,  the  range  output  will  be  in  feet.  Use 
this  option  only  for  a short-range  trajectories. 

Line  3 (Format  10F5.0):  The  first  field  contains  WO,  the  loaded 
(launch)  weight  of  the  vehicle,  including  fuel,  boosters,  etc.  The 
second  field  contains  WE,  the  vehicle  empty  weight,  without  usable  fuel, 
dropables,  etc.  Weights  must  be  in  pounds.  The  third  field  contains 
AREF,  the  aerodynamic  reference  area,  in  square  feet.  The  fourth  field 
contains  the  angle-of-attack  limit  in  degrees.  The  program  does  not  use 
negative  angles  of  attack,  and  a lower  limit  of  0 degrees  is  assumed. 

The  fifth  field  contains  AC,  normally  the  inlet  capture  area  in  square 
feet  per  sustainer  engine  of  unit  size.  AC  is  automatically  scaled  by 
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Figure  6.  Batch-Mode  Output  Format 
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the  input  data  value  TM.  AC  may  be  used  for  another  purpose  if  inlet 
capture  area  has  no  bearing  on  the  trajectory,  since  it  is  included  in 
the  common  block..  The  sixth  field  contains  TM,  the  thrust  multiplier  or 
sustainer  engine  scaler.  TM  scales  the  sustainer  thrust,  inlet  drag, 
airflow,  etc.  The  seventh  field  contains  GPU,  the  normal  load  factor 
used  for  pull-up  maneuvers  performed  in  M5 , in  g's.  The  last  three 
fields  contain  the  launch  parameters  XLH  (launch  altitude),  XLMN  (launch 
Mach  number),  and  XLFEE  (flight  path  angle  above  horizonal),  respectively. 
All  input  altitudes  are  in  thousands  of  feet  MSL,  and  all  input  angles 
are  in  degrees. 

Line  4 (FORMAT  11F5.0):  The  first  field  contains  DFEEC,  the  desired 
or  initial  climb  flight  path  angie  for  use  by  M2  and  M7.  The  second 
field  contains  CLMN,  the  desired  climb  Mach  number  for  use  in  M2.  If 
the  M-sequence  begins  2383...,  attainment  of  CLMN  in  the  first  occupation 
of  M3  will  trigger  transition  to  the  M8  climb.  The  third  field  contains 
BAC,  the  vehicle  body  angle  to  be  used  for  climbs  in  M8.  The  fourth 
field  contains  BAD,  the  body  angle  to  be  used  for  M8  dives.  The  fifth 
field  contains  DCH,  the  desired  cruise  altitude.  This  altitude  will  be 
sought  during  the  first  climb  or  dive  if  3 is  included  in  the  M-sequence. 
The  sixth  field  contains  CMN,  the  desired  cruise  Mach  number,  for  use  in 
M3.  The  seventh  field  contains  RC.  This  is  the  range  in  nautical  miles 
which,  when  attained,  causes  a transition  to  the  next  M-subroutine  in  the 
M-sequence.  Normally  RC  triggers  the  end  of  cruise.  When  the  M-sequence 
includes  ...310...,  RC  is  assumed  to  be  the  distance  from  launch  point 
to  target.  The  eighth  field  contains  DFEEDD,  the  desired  flight  path 
angle  for  use  in  M4  and  M9.  For  dives  and  landing  approaches,  positive 
values  should  be  used.  The  ninth  field  contains  DVMN,  the  Mach  number 
sought  in  M4  and  M9.  The  tenth  field  contains  DTH,  the  altitude  to  be 
sought  and  held  for  cruising  in  M5.  DTH  is  also  the  field  elevation  for 
M9  landings.  The  last  field  contains  DMN,  the  Mach  number  to  be  sought 
and  held  in  M5. 
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SUBROUTINE  FLYIT3  LISTING 


NWC  TM  3540 


HC  SUBROUTINE  FLYIT3  (MAIN  PROGRAM) 

-2«C  INPUT  INSTRUCTIONS  

3«C  A COMPLETE  DATA  DECK  CONSISTS  OF  3 CARDS  PER  TRAJECTORY 

~4:C  - FIRST  DATA  CARD  - TRAJECTORY  PHASE  SEQUENCE  IN  1011  FORMAT,-'  M~VAUUES 

5«C  M= I RUNWAY  TAKE-OFF  AT  MAX  POWER  FROM  FIELD  ELEVATION  XLH 

-<5«C  - M=2  AIR  LAUNCH  AT  XLH  ALTITUDE,  XLM  MACH,  XLFEE  FLIGHT  PATH' ANGLE,  AND 

7«C  CLIMB  AT  CLMN  MACH  DFEEC  FLIGHT  PATH  ANGLE  AND  CONTINUOUSLY  SEDUCE 

-8*C  — --FLIGHT  PATH  ANGLE  IF  NECESSARY  TO  MAINTAIN  SPEED.  

9 « C M=3  CRUISE  AT  DCH  ALTITUDE  AND  CMN  MACH  UNTIL  RC  RANGE  IS  REACHED 

-10iC OR  FUEL  IS  EXHAUSTED 

I 1 «C  M=4  DIVE  AT  -DFEEDD  FLIGHT  PATH  ANGLE  UNTIL  PULL-UP  IS  INITIATED 

1 2 *C  - - M=5  PULL-OUT  FROM  DIVE  AT  GPU  LOAD  FACTOR  AND  CRUISE  AT  DTK  ALTITUDE  AND 

I3*C  DMN  MACH 

1 4 «C  — M=6  FLY  AT  CONSTANT  VELOCITY  UNTIL  DCH  ALTITUDE  OR  RC  RANGE  IS  REACHED 

I 5«  C OR  IMPACT 

"1 6 * C — - M® 7 —BOOSTER  LAUNCH  WITH  FLIGHT  PATH  ANGLE  CF EEC  FIRST  ICO  FEET 
1 7 * C M=8  FLY  AT  BODY  ANGLE  3AC  OR  BAD  DEGREES  UNTIL  DCH  IS  REACHED 

I 8*C  - M=9  APPROACH  TO  LANDING  AT  GLIDE  SLOPE  DFEEDD,  RUNWAY  ELEVATION  DTH,  ~ 

1 9 « C APPROACH  SPEED  DMN  MACH 

. 20:C  — MslO  DIRECTS  THE  VEHICLE  TO  FLY  ALONG  THE  RADAR  HORIZON  WHEN  THE 

21 «C  HORIZON  IS  INTERCEPTED  IN  M=3  CRUISE.  TARGET  ANTENNA  IS  ASSUMED 

"22>C TO "BE  AT  100  "FEET  ABOVE  SEA  LEVEL.  RANGE  To  TARGET  MUST  3E  INPUT' 

23 tc  AS  RC  (NMI). 

24tC  — IF 'USED  AT  ALL,  M=l  AND  M=7  MUST  BE  FIRST  IN  THE  SEQUENCE. 

25  »C  IF  USED  AT  ALL,  M=9  WILL  BE  LAST  IN  THE  SEQUENCE.  _ 

26 :C  " — NO~M  VALUE  MAY  BE  USED  MORE  THAN  ONCE  IN  A TRAJECTORY  EXCEPT  2,3, 4,6,3. 

27 » C Ms4  MUST  BE  USED  IN  THE  SEQUENCE  45,  OR  BE  THE  LAST  IN  THE  SEQUENCE. 

"28«C M= 1 0 IF'USED " AT  ALL  MUST  BE  LAST  IN  THE  SEQUENCE 

29«C  REMAINDER  OF  FIRST  CARD  IN  FORMAT  815.  -- 

30 i C BOOSTER  MULTIPLIER 

31 «C  ATMOSPHERE  - 0 FOR  STD,  1 FOR  HOT  DAY,  2 FOR  COLD  DAY 

32, c AERODYNAMICS  - SUBROUTINE  NUMBER  FROM  AERO  LIBRARY 

33tC  BOOST  PROPULSION  - SUBROUTINE  NUMBER  FROM  BOOSTER  L I Bn  ARY  - C IF  NONE 

34 tC SUSTAIN  PROPULSION  - SUBROUTINE  NUMBER  FROM  SUSTAINER  LIBRARY  - 0 IF 

35 i C NO  SUSTAINER 

36 »C  - -"OUTPUT  OPTIONS  - 0 GIVES  BRIEF  FORMAT  (REPEATS  DATA  CARD  ENTRIES) 

37 «C  1 GIVES  FULL  PRINTER  LIST  OF  INPUT  DATA 

38, c GIVES  option  1 PLUS  A PLOT  OF  THE  TRAJECTORY 

39*C  3 INDICATES  DEMAND  TERMINAL  OUTPJT  FORMAT  (FULL  LIST 

-40»C "4  "INDICATES  SHORT  FORM  DEMAND  TERMINAL  OUTPUT 

41 :C  5 INDICATES  OPTION  4 PLUS  TRAJECTORY  PLOT 

42»C  CLOSER-LOOK  OPTION  - IF  FINER  PRINT  RESOLUTION  DESIRED,  PUNCH  1 IN  45 

43 «C  RANGE  OUTPUT  WILL  BE  IN  FEET  IF  I IS  PUNCHED  IN  50 

44:C  SECOND  DATA  CARD  - MISSILE  LAUNCH  WEIGHT  INCLUDING  BOOSTER,  MISSILE  EMPTY" 

45 «C  WEIGHT  INCLUDING  UNUSEABLE  FUEL,  MISSILE  REFERENCE  AREA  IN  SQ  FT, 

46 1 C " —ANGLE  OF  ATTACK  LIMIT  IN  DEGREES,  ENGINE  INLET  CAPTURE  AREA  IN  SQ  FT,- 
47«C  SUSTAINER  THRUST  MULTIPLIER,  LOAD  FACTOR  TO  BE  USED  IN  ANY  PULL-UP 

48 «C  MANEUVER,  LAUNCH  ALTITUDE,  LAUNCH  MACH  NUMBER, LAUNCH  FLIGHT  PATH  ANGLE 

49 i C - ALL  IN  FORMAT  I2F5.0. 

50  »C  THIRD  DATA  CARD  - DESIRED  CLIMB  FLIGHT  PATH  ANGLE,  CLIMB  MACH  NUMBER,  CLIMB 
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BODY  ANGLE,  DIVE  BODY  ANGLE,  CRUISE  ALTITUDE,  CRUISE  MACH  NUMBER,  CRUISE 
RANGE,  DIVE  ANGLE,  DIVE  MACH  NUMBER , DASH  ALTITUDE,  DASH  MACH  NUMBER 

- ALL  IN  FORMAT  I2E5.0.  _ _ __  

INPUT  DATA  UNITS  IN  DEGREES,  THOUSANDS  Or'  FEET,  SQUARE  FEET,  MACH  NUMBER, 

_ NAUTICAL  MILES,  G'S,  POUNDS  

MAIN  PROGRAM 

DIMENSION  M (10) , PP ( 2000 ) , RS ( 900 ) , HS ( 900 ) , RSFTt  400) , HSLD( 400 ) 

COMMON  W, WO, WE, R(),AREF,CLAR, ALPHA, ALPHA! , ALPL, V , VR , VH, Q.THH, FEED,  

1 DCH ,DT,H, DFEEC, 04, AMACH, CMACH , R , RC, GPU , JCT, DFEEDD,DFEADD, LCT , TI T, 

2NTD, TIME, PEER, ALPHR.D.XLH, BAG, BAD, ITD.DT3, I KT.NB, HS.HS, SSFT.CL, CD, 

3HSLD, MCT, NCT , TM , AC, CMN , CLMN , DVMN , DMN , PT , PRI , TT, ASPEED, AMACH I , FI SP, 

4THRS,  THRB  ,T , CDO  ,CDI  , CD  A , CDAOFF,  MM,  M,  PP , MGO , DTH , XLFEE,  XL.MN 

5,P, A, FISPB.FI SPS,NBSTR,NSUST,NAERO, BA 

THE  FOLLOWING  FUNCTION  IS  A LINEAR  INTERPOLATION  

TERPI ( X , A , B ) = B ♦ (A  - B)  * X 

THE  FOLLOWING  FUNCTION  INTERPOLATES  ALONG  AN  EXPONENTIAL  DECAY  (NON-LINEAR). 

A IS  VALUE  OF  HIGHER  ORDINATE,  B IS  THE  ABCISSA  OF  INTEREST  BETWEEN  C AND  E, 

_ C IS  THE  LOWER  ABCISSA  VALUE,  D IS  THE  . LOWER  J)RD  I NATE..  VALUE,  AND  E IS  THE 

HIGHER  ABCISSA  VALUE. 

TERP2 ( A ,3 ,C , D, E )=A*EXP ( ( B-C) *LOG ( D/A ) / ( E-C ) ) 

THE  FOLLOWING  FUNCTION  IS  A LINEAR  INTERPOLATION 

TERP3CXI  , X2, X , Y 1 , Y2.)_=Y I + ( Y2~Y_1  J *( X-X I )/( X2-XI.1 

WRITE (6, 20) 

_20  FOR MAT ('  ENTER  DATA.') 

READ  INPUT  DATA  CARDS 

■ I READC5,  2,END=9999)M(  I ) , MC  2)  , M( 3 ) , M(  4 ) , M ( 5) . M<  6 ) , M<7 .).,«(  3).,M.C2J., 

IM( 10) , NB, N ATM, N AERO, NBSTR , NSUST , NOUTPT, NPRT.NFEET 

_ 2 FORMATt  1011  ,815)  _ ..  

READ(5,3)WO,WE, AREF, ALPL, AC, TM, GPU, XLH.XLMN, XLFEE 

3 FORMAT! I2F5.0)  

REA DC  5 ,3 )DF EEC, CLMN , BAC, BAD.DCH, CMN,RC, DFEEDD, DVMN, DTH, DMN 

IF(NOUTPT.NE.O. AND. NOUTPT. LT.4)G0  TO  ICO  

WRITE! 6, 1 2) (M( I ) , 1= I , I 0 ) , N B, N ATM, N AERO, NBSTR, NSUST, NOUTPT, NPRT, 

I NFEET  

12  FOR MAT <y I', I 011  ,815) 

___WRITE(6,I3)W0, WE, AREF, ALPL, AC, TM, GPU, XLH.XLMN, XLFEE  _ 

13  FORMAT ('  ',F5.0,F5.0,F5.2,F5.0,F5.2,F5.2,F5.I ,Fb.2,F5.2,F5. I ) 

WRITE (6, I 4) DFEEC, CLMN, BAC, BAD, DCH, CMN, RC, DFEEDD, DVMN, DTH, DMN 

14  FORMAT ( F6 . I .F5.2.F5. I ,F5. I ,F5. I ,F5.2,F5.0,F5. I .F5.2.F5. I ,F5.2) 


DEFINE  SPECIAL  FUNCTIONS  AND  INITIAL  CONDITIONS 

100  XLH=XLH*I 000. 

DTH=DTH*I 000. 

DCH=DCH*I 000. 

NATM=NATM+ I 

MM=>  I 

ITD=0 

ECT=0 

JCT=0 
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-FLYIT3  CONTINUED 


01  « IKT=0  — 

02i  KCT=0 

-03* KP=  I 

04*  RS( I )=0.0 

105* RSFT(  i ) =0 .0 

06*  HS(I)=XLH 

,07*  — . — HSLD<  I )=XLH 

08*  K5=0 

; 09  * LCT=0 

110*  MCT=0 

I I I * NCT=0 

112*  NTD=0 

M3*  — - EWB=0 .0 

114*  KK=0 

H5,  -TIME=0.0 

M6*  DTB  =0.025 

117*  - TIM=0.0 

119*  R=0.0 

119*  - CDA=0.0 

120*  ASPEED= I 00. 

21  * fRIC=0.5 

122*  CMACH=0. 25 

(23,  H=XLH  

124*  FEED=XLFEE 

125* 1 F ( M ( MM ) . EQ . I )FEED=0.0 

126*  DFEEDD=-DFEEDD 

t-27  * DFEADD=DFEEDD 

128*  FEER=FEED/57. 2958 

129*  AMACH=XLMN  

130*  IF(M(MM).EQ.I)AMACH=0.0 

131  * - DO  8 MI  = I ,10  — 

132*  I F ( M ( M I + I ) . EQ . 1 )M(MI+I )=I0 

1 33 * 8 CONTINUE 

134*  Q3=5I8.688-0.00356616*XLH 

135*  V=AMACH*49.0I*SQRT(Q3>  - - 

136*  Q=0.00I*V**2.0 

137*  VR=V*COS( FEED/57. 3)  

138*  VH= V*S INC  FEED/5 7.3) 

1-39* AR I =0 . 0 

140*  AH  I =0.0 

(41*  VH I =0.0 — — 

142*  VRI =0.0 

143*  - I F ( M ( MM  ) -2 ) 34 , 33 , 33 

144*  33  ALPHA=6.0 

1-45* ALPHA  I =ALPHA 

146*  GO  TO  35 

147*  34  ALPHA=ALPL  

148*  35  ALPHR=ALPHA/57 . 3 

149*  AL=0. 9*W0  

150*  PR  I = I . 


& 
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FLYIT3  CONTINUED 


151 « 
I52i 

THR=AL/5.0  ’ - - ' 

THRB=0.0 

1 53 * 

1 54 « 

- - THRS=0. 0 * * ' --  --- 

FISPB=0.0 

1 55 1 

1 56* 

" FIS  PS =0.0 

THR()T=  1 .0 

1 57 » 
158* 

W=WO 

: 59* 
160* 

NLI NE=46 

NPRNT=4 

61  * 

; 62* 

CDX=0.  ~ * 

I F ( HOUT 3T . GT . 2 ) NPRNT=8 

163* 

164* 

' I F ( NPRT.GT. 0 ) NPRNT=NPRN'T/4  ” *" 

I F ( NO UTPT . GT . 2 ) NL I NE=50 

r 65  *’ 
166* 

IF! NOUTPT.LT. 1 ) GO  TO  I 01'  " " ' ' ' 

IFCNOUTPT.CT.O. AND. NOUTPT.LT. 4JCALL  0UTPT1 

167* 

168* 

101  GO  T0C102, 104,106), NATM 

102  WRITE (6, 103) 

169* 

1 70* 

103  FORMATS  STANDARD  ATMOSPHERE'/)  - 

GO  TO  30 

,71  « 
172* 

104  WRITE16, 105) 

105  FORMAT!'  HOT  DAY  ATMOSPHERE'/) 

f73* 

174* 

GO  TO  30 

106  WR ITE! 6 , 107) 

175* 

! 76* 

107  FORMAT!'  COLD  DAY  ATMOSPHERE'/) 

30  IF(NOUTPT.G£.3)WRITE(6, I 5) 

L7D  F 

C77* 

1 78* 

15  FORMAT!'  TIME  ALT  MACH  RANGE  WEIGHT-  THRUST  DRAG  ALPHA 

1 FPA  CDX  M'> 

'79* 

180* 

~7  R NO  UTPT . GE . 3 ) WR'  ITE!3 , 16 ) 

16  FORMAT!'  SEC.  FT.  NO  NM  LB  LB  LB  DEG 

DEG 

>81  * 
182* 

1 DEG  ')  " 

17  FOR MAT! F7. 1 ,F7.0,F5.2,F6. 1 ,F7. 1 ,F3. 1 ,F7. 1 ,F5. 1 ,F5. I ,F6. 1 , F5.2, 

13) 

183*  ‘ 

1 84* 

IF!N0UTPT.GE.3)G0  T036 

WRITE!6,3I ) 

1 85* 

> 86* 

J!  FORMAT!  I32H  -TIME VEHICLE  PERFORMANCE- 

1 MOSPHERE — -VEHICLE  ANCLES-  AERODYNAMICS -ENGINE 

--AT 

PERF 

87*“ 

88* 

' 20RMANCE-  M)  “ 

WRITE! 6, 32) 

89* 

' 90* 

32  FORMAT! I29H  SEC  MACH  HEIGHT  VH  AH  RANGE  VR  AR  WEIGHT  T A 

1 TT-2  0 BA  FP  ALPHA  LIFT  L/D  CD!)  CDI  DRAG  THRUST  RPM  T 

91  * 
192* 

2SFC  TTT7 ) 

36  DO  97  J=l , 19 

93* 

94* 

DO  91  “N= 1 ,NLINE 

DO  4 NN=I ,NPRNT 

r95* 

96* 

DO  94  1=1 ,5 

I F ( M ( MM ) . EQ . 0 ) MM=MM- 1 

97*- 

98* 

?XGO=M ! MM) 

MMI =MM 

99*  " 
!00« 

GO  T0!6,5,6,5,6,6^,'5,5V6,5),MG0 

5 IF! I.EQ. 1 )G0  TO  7 

A- 5 


30 


?0I 
!02 
>03 
>04 
>05 
>06 
>07 
>08 
>09 
>10 
>11 
>12 
>13 
>14 
>15 
>16 
>17 
2 18 
>19 
>20 
>21 
>22 
>23 
>24 
?25 
226 
>27 
>28 
>29 
>30 
>31 
232 
>33 
>34 
>35 

236 

237 

238 
>39 
>40 
>41 
242 
>43 
>44 
>45 

246 

247 
>48 
249 
>50 


GO  TO  I I 2 

6 I F ( I*NIJ.NE. I )G0  TO  1 12 

7 GO  T()(  109,110,11)  ) , NATM 

109  CALL  ATMSTD 

GO  TO  I I 2 

NO  CALL  ATMHOT 

GO  TO  I 1 2 _ 1 

111  CALL  ATMCLD 

J I 2 P=1 1 .91 46*S0*T 

“A=49.0I *SQHT(T) 

TT=T*( I .0+0.2*AMACH**2.0)  

TA=T-460. 

TRAT=TT/T  ; 

" PT=P*TRAT**3.5  ' 

_ 0=0 . 5*RO*V**2 .0  

DT=0. 25  ' ' 

.1121  GO  T0(  1 13,1  U,  I 15. 1 16, 1 17, 1 18, 1 19,1  20, 121 ,1  22), MGO 

113  CALL  Ml 

GO  TO  845  

114  CALL  M2 

GO  TO  845 

115  CALL  M3 

GO  TO  845  1 _ 

116  CALL  M4 

GO  TO  845 

117  CALL  M5 

J _G0  TO  845 ; 

118  CALL  M6 

GO  TO  845  

119  CALL  M7 
GO  TO  845 

120  CALL  M8 

GO  TO  845 

1 21  CALL  M9 

GO  TO  845  _ __  _ _ _ 

122  CALL  Ml  0 

845  IF( ALPHA. LT.O.O) ALPHA=0.0  _ 

I F( MM.NE.MM I )G0  TO  844] 

I F ( ALPHA. GT. ALPL ) ALPHA=ALPL 

ALPHA! “ALPHA 

842  I F C H ) 843 , 843 , I 23 

843  WHITE (6, 844) 

WRITE (6 , 846 ) 

~ WRITE! 6, 847 ) 

WRITE(6,848)  

WRITE(6 , 849  )*" 

844  FORMAT!/'  S _S  S SSSS  S SSSSS  SSSSS  _SSSSS  __S S !') 

1 !') 

846  FORMAT  ( ' SS  S ? S S S _ S S_S  __S_S S S !'  ) 

I') 
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-S-  --S S-  S -S S-S S — S 

S -SSSS—  SSSSS  - SSSSS  SSSSS  -SSSSS 


-FLYIT3-C0NTINUED- 


25 1 * ^ 847  FOR  MAT  C-7  SSS  SSSSS  -SSSS-  S S S_S S-SSSS- 

>52«  I * ) 

>531—848  FORMATS  S-S S S— S S-S -S S-S S — S 

>54«  I') 

>55*  849  FOR  MAT  ('  S S— S-  S -SSSS-  SSSSS-- SSSSS  - SSSSS— SSSSS 

>56*  \'/) 

>57*  - H=R+H/TAM ( F EER ) /6076 . 1 

>58*  RFT=R*6076.I 

>59* H=0 . 0 1 

>60*  KP=KP+I 

>61*  - RS(KP)=R 

>62*  HS(KP)=H 

>63*  GO  TO  92  — — 

>64*  8441  MCO=M(MM) 

>65* MM1 =MM 

>66*  GO  TO  1121 

>67 *C  LIFT  AND  DRAG  BLOCK 

>68*  123  GO  TO  ( 124, 125, 126, 127, 123, 129, 130, I 31 , 132, 133, 1331  , 1332, 

>69*  — I —1334,  1335), NAERO 

270*  124  CALL  AERO  I 

>71* GO  TO  134 

>72*  125  CALL  AER02 

>73* CO  TO  I 34 

>74*  126  CALL  AER03 

>75* GO- TO  134 

>76*  127  CALL  AER04 

>77*. — GO  TO-  1 34 

>78*  128  CALL  AER05 

>79*  GO  TO  134  

>80*  129  CALL  AER06 

>81*  - - GO  TO  134  

>82*  130  CALL  AER07 

>83* GO  TO-1  34- 

>84*  131  CALL  AER08 

>85*  - GO  TO  134  

>86*  132  CALL  AER09 

>87* GO  TO  134 

>88*  133  CALL  AERO  10 

>89  * GO  -TO  -134 

>90*  I 331  CALL  AERO  I I 

>91  * GO  TO  134 

>92*  1332  CALL  AERO  12 

>93*  GO  -TO -I  34 

>94*  1333  CALL  AERO  13 

>95* GO -TO-1 34 

>96*  1334  CALL  AERO  I 4 

>97*  — GO- TO  134  — 

>98*  1335  CALL  AER015 

>99*  134  IFCNBSTH.LT. 1 )G0  TO-145 

100*  GO  T0< I 35, 136, 1 37, 1 38, 1 39, 1 40, 1 41 , 1 42, 143, I 44, 1440, I 441 , 1 


1333, 


442, 1 443, 
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FLYIT3  CONTINUED 


301 

302 

303 

304 

305 
306 » 
307i 
308* 
309* 
310* 
31  I * 
312* 
313* 
314* 
315* 
316* 
317* 
318* 
319* 
320* 
321  * 
322* 
323* 
324* 
325* 
326* 
327* 
328* 
329* 
330* 
331* 
332* 
333* 
334* 
335* 
336* 
337* 
338* 
339* 

340  r 

341  * 
342* 
343* 
344* 
345* 
346* 
347* 
348* 
349* 
350* 


1 1 444 ) , NBSTR 

135  CALL  BSTHI 

GO  TO  145 

—136  CALL  BSTR2 

GO  TO  I 45 

- I 37  CALL  8STR3  — 

GO  TO  I 45 

-138  CALL  BSTR4  

GO  TO  145 

—339  CALL-BSTR5 

GO  TO  145 

140  CALL  BSTR6 — 

GO  TO  I 45 

"143  CALL  BSTR7  

GO  TO  145 

' M2  CALL  BSTR8 • ■ 

GO  TO  145 

143  CALL  BSTR9  

GO  TO  145 

144  CALL  BSTRIO  

GO  TO  145 

1440  CALL  BSTR11 

GO  TO  I 45 

1441  CALL  BSTRI2 

GO  TO  145 

3 442  CALL  BSTR13 r 

GO  TO  145 

T443  CALL  BSTR14  

GO  TO  145 

1444  CALL  BSTRI5  

145  IF( NSUST.LT. I ) GO  TO  1000 

GO  TO(l 46, 147, 148, 1 49, 1 50, 151 ,152,153,154,155,156,157,158,159,160, 
I I 61  , 1 62 , 1 63 , 1 64 , 1 65 ) , NSUST 

146  CALL  SUST1  

GO  TO  1000 

147  CALL  SUST2  ~ 

GO  TO  1000 

148  CALL  SUST3  1 

GO  TO  1000 

— I 49  CALL-SUST4 

GO  TO  1000 

-150  CALL  SUST5  

GO  TO  1000 

-151  CALL  SUST6 

GO  TO  1000 

152  CALL  SUST7  

GO  TO  1000 

153  CALL  SUST8  ‘ 

GO  TO  1000 

154  CALL  SUST9  
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GO  TO  1000 

155  CALL  SUSTIO  

GO  TO  1000 

- —1 56  CALL  SUSTII 

GO  TO  1000 

15-7  CALL  SUSTI2 — 

GO  TO  1000 

_.158  CALL  SUST13 

GO  TO  1000 

—159  CALL  SUSTI4 

GO  TO  1000 

- 160  CALL  SUSTI5 

GO  TO  1000 

- 1 61  CALL  SUSTI-6 

GO  TO  I 000 

-162  CALL  SUST1 7-  

GO  TO  I 000 

• 163  CALL  SUST18 

GO  TO  I 000 

164  CALL  SUST19- : 

GO  TO  1000 

-165  CALL  5UST20 

1000  CDO=CDO+CDAOFF 

CDX=CDO  

IFCXLFEE.GT. 0.0001 . AND. XLFEE. LT. I 0. )CDX=XLFEE 

— SFC=3600./(  THRS+THRB  )-*(THRB/F  ISPB+THRS/F  ISPS) 

FISP=3600./SFC 

THR=THRB+THRS 

: FLIGHT  EQUATIONS  OF  MOTION 

0=0 . 5*R0*  V **2 . 0 — 

CD-CDO+CDI 

W=W-DT*THR/F ISP/2.0 

D=AREF*CD*Q 

AL=AREF*CL*Q 

ALD=CL/CD 

ALPHR=ALPHA/57.3 

FR=THR*COS( FEER+ ALPHR ) - AL*S I H ( FEER ) -Q*COS  ( F EER ) 

- - -FH=THR*S I N ( F EER+ ALPHR ) + AL*COS ( F EER ) -D*S I NC  F EER ) -W 

IF( ITD.LT. I )G0  TO  1007 

FR=FR-N*FRIC 

1007  AR=32 . I 725*FR/N 

- - IF(ITD-I)I004, 1001, 1001  

1004  I F ( M C MM ) - 1 11002,1001 ,1005 

-1001  AH=0.0  

VH=0. 0 

AH1 =0.0 

GO  TO  1003 

1 005  IF(  MG0.NE.7  )G0  TO  1 002 

AA=32 . I 74*( THR-D-W*S INC  FEER ) ) /W 

AR=  AA*COS ( F EER ) 
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FLY IT3  CONTINUED 


101 

1C2 

103 

!04 

105 

lt>5 

107 

t08 

109 

110 
1 1 1 
112 

113 

114 

115 

116 

117 

118 

119 

120 
121 
122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 


I 002 
1003 


AH= AA*S IN  C F EER ) 

GO  TO — 1 003 

AH=32 . I 725*FH/N 

VR=VR+AR  1 *DT+0.  5*< AR-AR I)  *DT 

VH=VH+AHI*DT+C.5*(AH-AHI )*DT 
-R=R+(VR1 *DT+0.5*<VR-VR1 )*DT 1/6076. 
H=H+( VH1 *DT+0 ,5*( VH-VHI )*DT) 

ARI=AR 1 

AHI =AH 


-VR1=VR 

VH1 =VH 

FEER=ATAN ( VH/VR 1 

FEED=57 . 2958*FEER 

01 =VH**2.0+VR**2.0  

V=SQRT(QI ) 

" AMACHi  =AMACH 

AMACH=V/A 

W=W-DT*THR/F I SP/2VG  ' - 

GA=(THR-D)/W 

GL=AL/W — 

BA=ALPHA+FEED 

Q4=Q 

RTGT=RC-R 

I F ( M( MM+ 1 ) . NE . 1 0 . AMD . M ( MM+2 ) . ME . I 0 1 GO  TO  1006 
RHl)RH=.663 I 4*{RTCT-)2.28)**2 

“IF(H.GT.RHORH) MM=MM+ 1 7~ 

MGO=M( MM) 

‘ I F ( H . GT.'l  05 1 00 TOTGCTTO ~ 40 

IF( AMACH. GT.5.6 )G0  TO  40 

I F(H) 843, 843, 94 

TIME=TIME+DT 

CONTINUE  

RFT=R*6076.  I 

R0C=VH*60. 

I F( I . NE. 3. AND. I .ME. 5 )G0  TO  940 

KP=KP+ 1 - ‘ 

RSCKP)=R 

HS(  KP)=H  : 

940  I F( NBSTR.LT. I )G0  TO  92 

IFCTIME. GT. 10.) GO" TO  92  

RSFT ( N+ 1 )=R*6076. 1 

HSLD(N+1 )=H  

I F ( NOUTPT  .GT.  2 ) GO  TO*  93 


1006 


94 
- 4 


92 

98 

99 


93 

930 


WR I TE ( 6 , 99 ) TI  ME , AMACH,  H , VH . GL , R , VR , GA , 17,  TA , TT  , Q,  B A , FEED,  ALPHA , AL , 
ALD,CDO,CDI , D,THR, ASP EED, SF C, TIT, MGO 

FORMAT!  F7.1  , F5. 2,  F7 .0.F6.0,  F4 . 1 , F5. 0,  F6“.0,F4. 1 ,-F6.0,F5“.0,F5.T),F6.0- 
,F6. I ,F6. 1 ,F5. 1 ,F6.0,F4. 1 ,F4.2,F4.2,F6.0,F6.0,F6. I ,F5.2,F6.0, 12) 

GO  TO  932  • ---  --  

I F ( NFEET) 930,930, 93 1 

WR  I TE  ( 6 , 1 7 ) TI  ME , H,  AMACH,  R , 17  ,THR , D,  ALPHA , ALD , F EED , CDX , MGO  
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151  » GO  TO  932 

152 « 93 1 WHITE! 6, I 7)TI ME, H, AMACH, HFT.W ,THR,D, ALPHA, ALD, FEED, CDX , HGO  — 

153*  932  IFCH.LT.O.I )GO  To  40 

1541  - 91  -CONTINUE  — 

155*  95  WRITE(6, 19) 

156*  19  FORMAT!  I HI  ///////) — 

157*  I F ( NOUTPT.GT. 2 ) GO  TO  90 

*58*  - '96  WRITE (6, 3 I ) • 

159*  WRITE! 6, 32) 

160* CO- TO  97 

161  * 90  WR I TE ! 6 , 15) 

162*  97  CONTINUE  

163*  40  IF! NOUTPT. EQ. 2 . OH.NOUTPT. EQ. 5 )CALL  SHOMAT! PP, KP ,RS, HS ) 

164*  9999  CALL  EXIT 

165*  END 
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1*  - SUBROUTINE  OUTPTl-  _ - - _ 

2 * DIMENSION  M( 10), PP* 2000) ,RS( 900) ,HS*900) ,HSFT( 400) ,HSLD(400) 

3 *-  COMMON  N,WU,W E,RU,AR EE ,£LAR, ALPHA ,ALPHA1  , ALPL , V » VR ,-VH ,D ,IHR,£  EED  » 

4i  1 DCH ,DT, H,DFEEC, Q4, AMACH ,CMACH,R , RC, GPU , JCT, DFEEDD,DFEADD, LCT, TI T, 

5s  2NTD,TIME,FEER,ALPHH»D,XLH,BAC, BAD,ITD,DTB, IKT,NB,RS,HS,RSFT,CL, CD, 

6«  3HSLD, MCT, NCT , TM , AC , CMN , CLMN , D VMN , DMN , PT , PR  I , TT , AS P EED , A MACH  I , F I S P , 

7 s 4THRS , THRB ,T,CDO t CDI , CDA , CD AOFF , MM,  M , PP , MGO , DTH, XLFEE, XLMN 

3 s 5, P, A, FISPB, FISPS,NBSTR , NSUST, NAERO, BA 

9*-  WRITEC6-,6)W0 

10s  6 FORMAT ( 49H 1 GROSS  FUELLED  MISSILE  WEIGHT,  INCLUDING  BOOSTERS , F7 . I , 7 

lit  tH  POUNDS)-  — - 

1 2 s WR I TE ( 6 , 7 ) W E 

13*  7 F0RMATC39H  MISSILE  EMPTY-  WEIGHT,  WITHOUT  BOOSTERS  ,F7 -I  , 7H  POUNDS) 

14s  WRITE ( 6 , 8 ) AREF 

I 5 s. 8— FORMAT*  29H-DRAG-  AND— L-I F T- REF E RE NC a AREA,F6. 3,  1-2H  SQUARE-FEET) 

16*  WRITE (6, I 0)DFEEC 

17* 10  FOR  MAT*  39H-DESIRED- CLIMB- ANGLE  TO  CRUISE  ALTI TUDE  ,F5.-1-,SR  DEGREES) 

13$  WRITE* 6, 1 1 )CLMN 

19*. H -FORMAT*  4XH-DES I RED-CLIMB-  MACH  NUMBER,-  IF  APPLICABLE,F5.-2)- 

20$  WRITE (6, I 2) DCH, RC 

2 1 *- 1-2  FOR  MAT-02444— DESIRED— CR  U I SE-ALT  ITUDErFB,Qr9H-FE&T-*-TOTF6.0f  T3H -MILES 


22$  I RANGE)) 

23* WRITE*  6,1-3 ) G MN-  — 

24*  13  FORMAT*  27H  DESIRED  CRUISE  MACH  NUMBER, F5. 2) 

25* WR  LIE*  6 r 4-4-IDF  EEDD 

26$  14  FORMAT* 43H  DESIRED  DESCENT  ANGLE  FROM  CRUISE  ALTITUDE, F6. 1 , 9H  DEGR 

27  * 1 EES, ) 

28$  WR I TE ( 6 , 1 5 ) DV MN 

29-t-  -15  FORMAT*  28H  DESIRED -DESCENT- MACH  NUMBER, F6,3-> 

30*  WRITE (6, 1 6) DTH  , DMN 

3is  16  FORMAT ( 22H  DESIRED  DASH  ALTITUDE, F8. 0, 1 I H FEET  *MACH, F5.2 ,2H  )) 

32*  WRITE(6,36)TM 


33* 36  FOR MATI47H- SUSTA I NER  THRUST-MULTIPLIER  (X  SASIC-THRUST-)— =-,F5.-2 ) 

34s  WRITE(6,5)NB 

35*-  5 FORMAT*'  BOOSTER  MULTIPLIER- = ',11)  — 

36s  WRITE (6, 28) AC 

37  $ 28  FORMAT*  27H- INLET  C*M)K  I E-CUTTER-  AREA  =,F8.3, 12HX— TM,-SO»FT^) 

38$  WRITE (6, 38) GPU 


39  J 


40*  WRITE (6 , 25 ) ALPL 

41  $—  2S- FOR  MATT  28H-  LI  MTT -ON-ANGLE-  OF- ATTACK -=-,F  4,3  ,SH  DEGREES) 

42*  WRITE (6, I )BAC,BAD 

43*  . - l - FORMAT C' - BODY- ANGLES,-  IF-APPLICABLS  - TN.  CLIMB'  ,F5,J  ,' IN-DI V t'  , F 

44 s 15.1,'  (DEGREES)-') 

45t XFX  M4  MM)=X4-l-T,L7^20 

46*  17  WRITE (6, 1 8) 

47*  18  FORMAT*  69H  THIS- MISSILE  IS  SURFACE-LAUNCHED- BY-  GROUND-  ROLL- UNDER 


48*  1 ITS  OWN  POWER.////) 

49*  GO  TO  29 


50*  20  WRITE* 6,21  ) 
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51 $ 21  FORMAT ( 81H  THIS  MISSILE  IS  AIR-LAUNCHED  AND  SUBJECT  TO  THE  FOLLO. 

52 * LING  ADDITIONAL  CONDITIONS.)  _ . _ 

53  s NR  I TE ( 6 , 22 ) XLH 

_ 54  s 22  FOR  MAT  (Lori  LAUNCH. ALTITUDE.F8.Q..5H  FEET )_  _____  

55s  WRITE (6, 23 >XLMN 

56 s 23  FORMAT.Q9H  LAUNCH  MACH .NUMBER *F5. 3) 

57$  NR I TEC 6, 24)  XLFEE 

58$  24  FORMAT ( 25H  LAUNCH  FLIGHT  PATH  ANGLE, F5. I ,5H  DEG.//) 

59$  29  'NR  I TEC  6, 26)  CMC  I ) , 1=1  , 10) 

60$  26 . FORMAT  t'  .TRAJECTORY  SEQUENCE  - .M  Oil//.)  

61  $ NRITE (6 ,2 ) NBSTR , NB, NSUST, NAERO 

62$  2 FOR  MAT C'  OPTIONS  CALLED  _ BOOSTER' , 13 NUMBER  OF  BOOSTERS  = ' 

63$  1 ,12,'  SUSTAINER', 13,'  AERODYNAMICS' , 13///) 

64$  IFCN0UTPT.EQ.2) WRITE (6*27) 

65$  27  FORMAT C ' TRAJECTORY  PLOT  REQUESTED'///) 

66$  NRITE  C6 ,19) . 

67$  19  FORMAT C 1 HI ) 

68$  RETURN . - __ 

69  $ END 


B-3 


NWC  TM  3540 


Appendix  C 
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c 
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1 


NWC  TM  3540 


X • 

< in  • x 

T*  7" 

XXX  - 

»> 

'X  • LU  LU 

Q 10  ’U  LJ 

o x 

• XXX 

o > 

X -in  j 

•* 

X X < X 

— > 

X Z • . 1 

o 
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SUBROUTINE  SUSTI  . ...  

DIMENSION  MC 10) , PR ( 2000) , RS( 900 ) ,HS ( 900 > , RSFT< 400) ,HSLD(400) 

COMMON  W , WO , WE . HO , AREF , CL AH , ALPHA , ALPHA  I , ALPL , V , VR , VH , Q , THH . F EED , 

1 DCH, DT , H, DFEEC, Q4, AMACH.CMACH , R , RC.GPU , JCT, DFEEDD, QFE ADD, LCT.TIT, 

2NTD, TIME, PEER , ALP HR,  D,XLH,B AC, BAD, I TD.DT3, I KT,NB, RS ,HS , RSFT,CL, CD, 

3HSLD, MCT, NCT , TM , AC , C MN , CLMM , DVMN , DMN , PT , PH  I , TT , ASPEED , A MACH  I , F I SP , 

_ 4THRS,THRB,T,  CDO , CD  I ,CDA,CDA()FF,  MM,  M,  PP.MGO,  DTH,  XLFEE,  XLMN 
5,  P,  A,  FI SPB, FISPS,NBSTR, NSUST, NAERO, 3 A 
C CURTI SS-WRIGHT  ETE  3 MAY  78  NORMAL  SHOCK  INLET  3 PERCENT  DUCT  LOSS 
C NO  BYPASS  NO  BLEED  THROTTLE  ENGINE  TO  SUBCRITICAL 

PR I =0.97 

IF(W-WE)70I ,702,702 

.70)  THRS=0.  . 

GO  TO  .799 

. 37  FORMAT ('  C-W  ETE  3 MAY  78  NORMAL  SHOCK  INLET  3 PCNT  DUCT  LOSS') 

702  IF(AMACH- I .0)703,703,704 

._703  DELTA=PT/I4.696*PRI 

GO  TO  705 

704  DELTA=PRI*PT/1 4.696/C I .0+1 .1 6667*(AMACH**2. 0- I .0 ) ) **2.5/( 0.  1 66667+ 

1 0.83333/AMACH**2. 0) **3.5 

705  THETA=TT/51  3.7  . — - _ 

I F(  TT-550 .0 ) 706,706,707 

—.706  R P MM AX=  ( 2 6 960  . + TT*3 . 2 72 72 73 )/ 2 86 . 5 

GO  TO  708 

— ' 707  RPMMAX=C26581 .+TT+3.96I 538)/286.5  

708  MGO-M(MM) 

KTHHT=0 - 

GO  TO  (715,715,709,710,71 1 ,715, 709, 715, 710, 710), MGO 

709  DAI =CMN 

GO  TO  712 

710  IF ( I TD)  7 1 7,7 1 7, 720 

7)7  DM=DVMN 

GO  TO  712  : 

720  DM=0.0 

GO  TO  712 

711  DM=DMM 

712  ASPEED=( ASPEED+5. *( DM-AMACH ) ) -I  00 ,0*( AMACH-AMACH1 ) 

GO  TO  716 

.715  I F ( M ( MM+ 1 ) . EQ  . 5 ) GO  TO  717  

ASPEED=RPMMAX 

KTHRT= 1 

716  I F ( ASPEED.GT. RPMMAX ) ASPEED=RPMMAX 

IF( ITD.GT. 1 ) ASPEED=0.95*ASPEED  - 

IF< ASPEED.LT. 80.0) ASPEED=60.0*SQRT( TMETA)  . 

...  IF(AMACH.LT.0.4.0R.ASPEED.LT.80. 1 )G0  TO.  7169  

I F < KTHRT.LT. I )G0  TO  7169 

I F( AMACH. LT. I . )AOAC=l .+1.1 5625*( 1 .-AMACH) 

7161  WDOTAC=-36 . 696+ .09568+DELTA+0 .985662 4238+RS PEED-. 0039 63 40333* 

I HSPEED++2 .0 


b4 


1 RSPEED**2 .0 

WDOT=i/D()TAC*DELTA*TM/SQfiT(  THETA)  - 

AO=0.09I  32*rtDOT/(R()*AMACH*SQRT(T)  )/l  44.0 

CAR=A0/AC/TM  

IF(  CAR.LT.A0AOG0  TO  7169 

ASPEED=ASPEED*£ I .-.85*1 1 .-AOAC/CAR))-.0I 

RSPEED=ASPEED/SCRT£ THETA) 

- GO  TO  7161  

7169  RSPEED=ASPEED/SQRT( THETA) 

— I F(H.GT. 20000.0) GO  TO  718  

IF(AMACH-0. 9)71 70,71 70,71 71  - 

-7170  AF=99. 4+0. 004420*H+ < -27  I .56-0. 000528*H)*AMACH 

BF=0. 81 6085-0.00001 08 1 *H+( 0. 1 930622+0. 00C00669U I *H) *AMACK 

— -GO  TO  7182 

7171  AF=-50.05+O.0l I 8925*H+ ( -I  05. 5-0. 008775*H) *AMACH 

QF=0. 7808925-0. C0000905*H+(0. 232  I 65+0. 0000047 43*H)*AMACH  

GO  TO  7182 

- -71  8 IFUMACH-I. 5)7180,7180,  7181 

7180  AF=226. 34-0. 001  95 2*h'+£-422.  45+0. 00707267*H)  *AMACH 

BF=0. 5228045+0. 00000385 4*H+(0. 500093-0. 0000086 53*H ) *AMACH  

GO  TO  7182 

7181  AF=41 0.07-0. 01  1 1 1 3*h'+(-536.  2+0.01  276*H)* AMACH  

BF=0. 5226605+0. 0C000386 1 6+H+ (0. 5001 89-0. 0000C86582*H ) *A MACH 

-7182  FD()TC=AF+BF*(  ABS  ( RSPEED-54.  I 1 ))**2.0 

7 1 9 FDOTA«TM*FDOTC*DcLTA*SQRT ( THETA ) 

— — THRDEL=( 48. 89-5 . 549+ ( ABS ( AMACH-O. 9) ) **3 . 3547 ) *RSP EED-C 3563. -523 .! * 

1 (ABS (AMACH-O. 9 ) )**2.5663) 

1 F( AMACH . LT.0 .9 )THRDEL= (47.5+1. 54444444+AMACH ) *RSPEED-( 3563 .-523.1- 

l*(ABS(AMACH-0.9 ) )**2.5663) 

— 722-THR=TM*0ELTA*THRDEL 

I F ( THR . LT . 0 . 0 ) THR  =+ I .0 

FISP=3600.-*THR/FDoTA 

I F ( F I SP . LT. 200 . CC ) F I S P=200 . CO 

TI TT2=0 .00003* ( RSPEED+ 1585.04) **2.0-80.629 

TIT=TITT2*rr 

1 F ( AMACH.  GT.  1 . )CDA=(  1 .517*  (LOG  (AMACH)  )**0.4)*(  1. 0-CAR  )-*AC*TM/AREF — 

I F(CDA.LT.0. 000001 )CDA=0.0 

CD()=CD()+CDA  — 

XLFEE=CDA 

THRS-THR  — — 

FISPS=FISP 

— 799- CDAOFF=  ( PT-+P ) *A C*66 . 0/SORT (-TT ) / ( 1 .0+0. 2*AMACH**2. 0)  **3.  0*S0RT(  5.0*_ 
1 ( (0.5*( PT+P )/P) **( . 286 )-l .0) )*(V-A*SQRT(5.0*( (0. 5*( PT+P ) /P) **( . 266 

2 )- 1 .0)  ))/32.  I 74*TM/0/AREF 

IF(THRS.GT.  0.01  )CDA()FF=0.0 

C. END  OF  SJSTAINER  DECK — 1 

RETURN 

END 
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SUBROUTINE  SUST2 

DIMENSION  M ( I 0 ) , PP ( 

COMMON  W,W(),WE,RO,AR 
“ 1 DCH , DT, H , DF  EEC , 04 , A !■' 
2NTD , T I ME , FEER , ALPHR , 
3HSLD , MCI,  NCT, TM , A C , C 
4THRS, THRB ,T, CCO  , CDI , 
5,P, A,FISPB,FIS?S 
THE  FOLLOWING  FUNCTION 
TER  P 1 ( X , A, B ) = B + ( 
THE  FOLLOWING  FUNCTION 


2000) ,RSC 900) ,HS( 900) ,RSFT( 4C0) ,HSLD(400) 
EF.CLAR, ALPHA, ALPHA  I , ALPL, V , VR, VH ,0 , THR , FE 
ACH,C.MACH , R , HC,GPU,  JCT,  DFEEDD,DFEADD,LCT,  T 
D.XLH, EAC , 3 AD , I TD , DT3, I KT , MB , RS , HS, RSFT, GL 
MN , CLMN , DVVN , DMN , PT, PR  I , TT, AS  PEED , A MACH  1 . F 
CDA , CDAOFF , MM , M , P P , MGO , DTH , XLFEE, XLMN 


IT. 

rn 

t '-lj  * 

ISP, 


IS  A LINEAR  INTERPOLATION 

A — B)~  * X 

INTERPOLATES  ALONG  AN  EXPONENTIAL  DECAY  (NON-LINE A 


C A VALUE  OF  HIGHER  ORDINATE,  B _IS_THE  *BCISSA  OF  INTEREST  BETWEEN  C AND  £, 

C C_IS  THE  LOWER !_ABCISSA  VALUE,_  D_IS  THE  LOWER  ORDINATE  VALUE.  AND  E IS_T!i£ 

C HIGHER  ABCISSA  VALUE. 

TERP2(A,3,C,D,E)=A*EXPrtB-C)*LOC(P/A)/(E-C) ) 

C PERFORMANCE  OF  J69-T406  TURBOJET  ENGINE  WITH  FULL  AFTERBURNER,  99  PCT 
C~ INLET  RECOVERY  AND  6 PCT  INLET  DISTORTION.  STD  DAY  ONLY.  4 DEC  74. 
IF(W-WE)70! ,702,702 

701~THRS=0.  

GO  TO  799 

37  ‘ FOR  MAT  ( >~PERF()R  MXNC E"T)F“J59-T40‘6~TOR Bo  J ET  ENGT  SET) ITH  TUCT  AFTER  B URN 

IRNER') 

102  MGO=M( MM) 

GO  TO  (715, 708, 709, 710, 711, 715, 715, 715, 710, 710), MGO 
708 "DM=CLMN  “ 

GO  TO  724 
709  ~DM=CMN 

GO  TO  724 

“7r0_ 

717 

720 

711 
~724 


715 

714 


IF(TTD)  7 17.717772CF 

DM=DVMN 
GO  TO  724 
DM=0.0 
' GO  TO  724 
DM=DMN 

ASPEED=(ASPEED+(0.0!0*(DM-AMACH)-0.050*(AMACH-AMACH! )/DT)*(H+ 

1 1 0000.0 ) )*0.5+0.5*ASPSED 

I F < AS PEED. GT. 100.0) ASPEED=I 00 .0 

IF(ASPEED.LT.90.0)ASPEED=90.0 

_G(rT()  714 

ASPEED= 100.0 


RSPEED=ASPEED*SCRT(  51 87688/TT)  “7 

AAA=(ASPEED*221  .5-1 9936. )/22 I 4. 

THR()T=ASPEED/ 1 00'. 

713  IF(H.LT. 20000. 0)G()  TO  719 

THR  40=  221  .4  + 652.1  *AMACH+I 95.2*AMACH**2-1  7702*AMACH**3 

THR40T=2224.8-31  99.9* AMACH+2 429. 5*AMACH**2. -4 78.94* A.MACH**3. 
SFC 40=2. 0284+0.  00601  *AMACH-0.0253*AMACH**2+0.04085*AMACH**3 
SFC40T= I . 746+0. 9303*A, MACH-O. 6792*AMACH**2.+0. 1 339*AMACH**3. 

IF( ASPEED.LT. ICO. )THR40=TERP1 ( AAA.THR40 ,THR 40T)  

IFUSPEED.LT.  100.  )SFC40=TERPI  ( AAA,SFC40,SFC40T) 
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SUST2  CONTINUED 


51  « 

52 
-53 

54 

55 

56 
-57 

58 

59 

60 
61 
62 

-63 

64 

-65 

66 

67 

68 
-69 

70 

-71 

72 

73 

74 
-75 

76 

-77 

78 

79 

80 
-81 

82 

83 

84 

85 

86 
-87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 


— 719  -THR20»I 443. 3-295. 4*AMACH+ I 327. 2+AMACH+* 2-272. 93*AMACH**3 

THR20T=  I 335.0-333 .2+AMACH+9 55 . 5*AMACH**2 .0- 1 53. 94*AMACH**3. 

SFC20=2 .06+0 . 1 905*AMACH— 0 . 2 I 593* A MACH** 2+0. I074*AMACH**3 

SFC20T=2. 186+0.2! 488+AWACH-O.  ! 692*AMACH**2.+0.07522*A.WACH**3. 

- -IF( ASPEED.LT. 100. )THH20=TE3P1  ( AAA, 7HR20,THR20T>  

IF( ASPEED.LT. I 00. >SFC20=TEHPI C AAA,SFC20,SFC20T) 

1 F < H . LT . 2 0000 . 0 ) GO  TO  718  

THR=TERP2 ( THR20 , H ,20000 . 0 , THR 40 , 40000 . 0 ) 

SFC=TERP2(SFC20, H, 20000. 0.SFC40, 40000.0)  

GO  TO  723 

718  THRSL  = 27 1 0. 4-68  I . 6*AM ACH+ 1 994 . 6* AMACH*+2— 33 1 . 7*AMACH**3  - - — 

THRSLT=2862. 1-3257.6*AMACH+4426.*AMACH**2.-I297.47*AMACH**3. 

SFCSL=2.079+0.50I06*AMACH-0.4442*AMACH**2+O.I  71  52 *A MACH** 3 

SFCSLT=  2.055+1. 5328+AMACH- 1 .5281 *AMACH**2. +0. 5237*AMACH **3 . 

IF( ASPEED.LT. 100. )THRSL=TERPI (AAA.THRSL.THRSLT)  

I F t ASPEED . LT. 1 00 . )SFCSL=TERP I ( AAA , SFCSL , SFCSLT) 

— - THR =TERP2(THRSL,H,0, THR 20, 20000.0)  

SFC=TEHP2 ( SFCSL, H.0.0.SFC20, 20000.0) 

-723  PRDUCT=0. 9898 

IF(AMACH.GT.  I .0  )PRDUCT=  (0.989  3*  ( 2 . 4*a,MACH**2. 0/  (0 ,4*AMACH**2 . 0+2 . 0 

- I ) )**3.5/( I .0+1 . I 6666667 *( AMACH**2 .0- 1 .0 ) ) **2. 5) / ( ! .0-0.075*1 AMACH- 
21 .0)**! .35) 

THR=THR*( I . 0-( I . 40+0 . 3903*( ABS ( I ! 2. 867-RSPEED) ) **4.0725  /I000000.0 

I )★( 1 . O-PRDUCT) ) 

DELTA=PT*PRDUCT/) 4.696 

THETA=SORT( TT/5 18.7) 

FISPS=3600. 0/SFC — 

THRS=THR*TM 

- '-VDOTAC=32. 1 -0.8382*(A3S(RS?EED*0.22l5-2  4.5)  )**!  .325 
TIT=TT*( I . 0+0. 0335* ( ABS (RSPEED^O. 22! 5-10.0) )**! .739)+350.0 

WDOT=WDOTAC*DELTA+TM/SQRT( THETA)  


A0=0. 09  ! 32*-'IDOT/(  R0-* AMACH*SCRT(  T ) ) / 1 44 . 0 

CAR=AO/AC/TM  - --  - 

CDA=0.0 

- IF(AMACH.GT. I .0)CDA=( 1.51 7* (LOG( AMACH ) ) *+0.4 ) *( I . O-CAR )*AC/AREF 
I F ( CD A. LT . 0 . 0 ) CDA=-0 . 00  I 

1 F(CDA)  725, 716,-716 


725  ASPEED=0. 999+ASPEED 
- GO  TO  714 
716  CDO=CDO+CDA 

799  CDAOFF= ( PT+P ) *AC+66 .0/SQRT( TT)/< I .0+0. 2* AMACH++2 
I ( (0. 5*( PT+P )/P) ★*< . 286 ) -I . 0 ) ) *( V-A*SORT ( 5.0*( (0, 

2)-! .0)  ))/32.l74*TM/G/AREF  

IFITHRS.GT.O.OI )CDAOFF=0.0 

END  OF  SUSTAINER  PERFORMANCE  DECK  - 

RETURN 


0) *+3. O+SQRT ( 5.0* 
5*(PT+P)/P)**( .286 


i 


G-5 


67 


NWC  TM  3540 


-SUBROUTINE  SUST3  

DIMENSION  M ( 10)  ,PP(2000) ,RS(900) ,HS( 900 ) t RSFT! 400 ) ,f!SLD(  400 ) 

— COMMON  W, WO, WE, RO, AHEF , CLAN  , ALPHA , ALPHA  I , ALPL,  V , MR,  VH , 0 ,THR , FEED, 

I DCH , DT, H, DFEEC, Q4, AMACH, CMACH , R , RC,  CPU , JCT,  DFEEDD , DFEADD,  LC'i  .TIT, 

-2NTD, TIME, FEES, ALPHR,D,XLH, BAG, BAD, ITD, DTB , I KT, NB, HS ,HS, RSFT, CL, CD,  — - 
3HSLD, MCT, NCT, TM, AC, CMN , CLMN ,DVMN , DMN , PT, PH  I , TT, ASPEED , A MACH  I , F I SP , 

4THRS, THHB ,T,CDO,CDI , CDA, CDAOFF, MM, M, PP, MGO, DTH, XLFEE, XLMN  

5,P,A,FISPB,FISPS, NBSTR, NSUST, NAERO, 5 A 
CONSTANT  THRUST  ROCKET  • MOTOR  WITH  232.9  ISP,  -19622  IMPULSE,- AND  380 -POON-S 

— IO*C  THRUST.  THRUST  MULTIPLIER  MAY  BE  USED  TO  CHANGE  THRUST  AND  CORRESPONDING 

— me --BURN  TIME  WHILE  MAINTAINING  TOTAL  IMPULSE. 

12*  702  THRS=0.0 

— 13* IF(  MGO.  EQ . 3 .OR.  MGO.  EO.  8 ) JCT=JCT+  I 

14*  IF< JCT.EQ.l )TIGN=TIME 

— 15*  IF<  JCT.LT.  I )G0  TO  799  

16*  THRS=380. *TM 

— I 7 * IF( (TIME-TIGN) .GT. ( 19622. /THRS) )THRS=0.0  

18*  FISPS=232 .9 

— 19*  -~799  CDAOFF=0.01  /AREF/AMACH- 

20*  IFI THRS. GT. I . ) CDAOFF =0.0 

— 21*  CDO=CDO+.02*TM/AREF  

22*  RETURN 

— 23* END 


1 , 

2* 

3, 

4* 

— 5* 

6* 

?«  

8* 

9*C  — 
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I 

— 2 
3 

4 

5 

6 

7 

8 

9 

— 10 
I 1 

— 12 
13 
— 14 
15 

— 16 
17 

— 18 
19 

- 20 
21 

— 22 

23 

24 

25 

26 
27 

— 28 

29 

30 

31 
-32 

33 

—34 

35 

36 

37 

38 

39 

- 40 
41 

- 42 
43 

- 44 
45 

— 46 

47 

48 

49 

50 


SUBROUTINE 
DIMENSION 


SUST4 

M< 10) , PP<  2000 ) ,HS( 900) , HS C 900 ) , RSFT( 400) ,HSLD( 400) 


COMMON  W , WO,  WE , so , ARE? , CLAR  , ALPHA,  ALPHA  ! , ALPL , V , VI! , VH 


I h!R , r El 


^ , 


701 


1 DCH , DT , H , D?  EEC , 04 , AM  ACH , CMACH , R , SC , GPU , JC f , DF EEDD , DEE A DD , LOT .TIT,- 

2NTD, TIME, FEES, ALPHR , D,XLH, 5AC , SAD, ITD, DTB, I KT , NB, RS ,HS, RSFT, CL, CD, 

3HSLD, MCT, NCT, TM , AC.CMN, CLMM , DVMN , DMN.PT, PR  I , IT, AS PEED, A MACH  I , FI SP , 

4THRS,  THRB,  T,  CDO  , CDI , CDA,  CDAOFF , M.M,.M,  PP,  MGO,  DTH,  XLFEE,  XLMN 

— 5,P,A,FISPB,FISPS, N8STR, NSUST, NAERO, BA  - — 

C THE  FOLLOWING  FUNCTION  IS  A LINEAR  INTERPOLATION 

— TER  PI  ( X , A , 3 ) = B + ( A—  B)  * X 

C TELEDYNE  CAE  J402-CA-400  HARPOON  ENGINE  21  JAN  75  GRB 

37  FORMAT!-'  J402-CA-400  MISSILE  TURBOJET  ENGINE  (HARPOON)')  

IF(W-WE)70I ,702,702 

THRS=0.  

GO  TO  799 

702  IFCAMACH-I .0)703,703,704  - 

703  DELTA=PT/I 4.696*PRI 

GO  TO  705  — 

704  DELTA=PSI*PT/1 4 . 696/ ( I .0+1 
I 0.8 3333/AM ACH**2.0)**3. 5 

705  THETA=TT/5 19,7 

713  RPMMAX=  I 00. 

GO  TO  <715,715,709,710,71 1 

709  D.M=CMN  

GO  TO  712 

710  I FC I TD ) 7 1 7,717,720 

717  DM=DVMN 

GO  TO  712 

720  DM=0. 0 

— GO  TO  712 

7 I I DM=DMN 

I F ( NTD . GT . 0 ) DM=0 . 0 

712  ASPEED=( ASPEED+ I 0.*(DM- 

GO  TO  716 

ASPEED=RP,MMAX 


I 6667* (AMACH**2. 0-1 .0) )**2.5/(0. 1 66667+ 


,715, 709, 715, 710, 710), MGO 


AM ACH) ) -500. 0* ( AM ACH-AMACH 1 ) 


7 l 5 
716 


IF( ASPEED.GT. RPMUAX ) ASP EED=RPMMAX 
I F ( ITD.GT.l ) ASPEED=0.95*ASPEED 
— I F ( AS PEED. LT, 80 .0 ) ASPEED=30.0  - 
RSP EED=ASPEED/SQRT( T1 1 ETA) 

— CN=408.*RSPSED  


706 

718 


730 


IF( RSPEED-I 00. ) 7 I 8, 7 I 3, 706 

RSPEED= I 00.0  

CN=408.*R SPEED 

ASPEED=RSPEED*SORT (THETA)  

I F ( AiMACH-0, 2)730,730,731 


CTHR I = • 
1 *CN**3 

CTHR2= 

I *CN*+3 
FDOTCl = 


.71 79E-4  — .9544264E- 
.81 70E-4  -. 7754703E- 


1 *CN+  .485031 75 E-5*CM**2— .51 85561 E- 10  — 
1 *CN+  . 3S823050E-5*CN**2  — . 3975  332E- I 0 


.3193E-3  -. 2292694E-I  *CN+  . 82364066E-6*CN**2+  .5I666I7E-I1 
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SJJ ST 4__CU  NT  IN  HE  D. 


51* l*CN**3  

52*  FDOTC2=  .2822E-3+  . 4430061 E-l *CN  26297909E-5*CN**2+  .4934756E-I0 

_ 53  *_  _I*CN**3  _ 

54*  X XX=AMACH/0 .2 

55*  _ GO  TO  721  

56*  731  IFUMACH-O. 4)732”,  732,733 

57*  732  CTHR 1 = .8I70E-4  7754703E-I *CN+  . 38823050E-5*CN**2  3975332E- I 0 _ 

58*  1 *CN  **3 

59* FDOTC I = 2622E-3+  . 44 3006 1E-I*CN  26297909 E-5*CN**2+  .493475 6 E-iQ 

CO*  |*CN**3 

6I_*_ CTHR2=  I757E-3_t.  I 035853E+0*CN+  .51  447335E-5*CIi**2  -,  55569_64E-  I 0 

62*  l*CN**3 

63 *_ FD0TC2=  .3668E-3  -.  9836357E-2*CN*  .1  IJ48998E-6*CN**2+.  _ .J_445J_9 4E-  1 0 

64  * I *cn **3 

65* XXX=(  AMACH-O.  21/0  ._2 

66*  GO  TO  721 

67*  J733  IFUMACH-0. 6)734,  .734, 735  _ 

68*  734  CTHR! = . I 757E-3  -. 1 035853E+0*CM+  . 51 447335E-5*CN**2  -. 5556964E- 1 0 

69  * 1 *CN**3 • 

70*  FDOTC 1 = . 3668E-3  — . 9836357E-2*CN+  . I I 748998E-6*CN**2+  .I445I94E-10 

71  * 1 *CN**3 

72*  CTHR2=  . 2357E-3  -. I 33 1 252E+0*CM+  . 66680686E-5*CN**2  -. 7552624E- I 0 

J 3* I *CN**3 

74*  FD0TC2=  . 525SE-3  -. I 234 748E+0*CN+  . 60270536E-5*CN**2  -. 6236748E-1 0 

75* 1 *CN**3  

76*  ’ XXX =( AMACH-O. 4 ) /O. 2 

77* GO  TO  721  

78*  735  I F( AMACH-O. 8 >736, 736, 737 

79*  736  CTHR!  = _ .2357E-3  -. I 33 1 252E+0*CN+  . 66680686E-5*CN**2  -. 7552624E-I 0 

80*  1 *CN**3 

81  «_ FD0TC1 = .5258E-3_-. I234  748E+0*CM+  ,60270536E-5*CN**2_r.6236  748E-l 0_ 

82*  | *CN**3 

83* CTHR2= .2324E-3  jtJ.3  1 3627E+0*CN+_  .66043l62E-5*CN**2_r,J75267_26t-J0 

84*  I *CN**3 

85*  __  FD0TC2=  .4769E-3  r.J739247E+0*CN+  . 87405830E-5*C?J**2_-.  9883  I 8OE1J  0 

86*  1 *CN**3 

87* XXX  = ( AMACH-O. 6)  /0.2  

88*  GO  TO  721 

89* 737  IFUMACH-I  .0)738,738,  739  

90*  738  CTHR I = .2324E-3  -.131 3627E+0*CN+  . 66043 1 62 E-5*CN**2  -. 7526726E-I 0 

9IJ l*CN**3 _ 

92*  FDOTC I = . 4769E-3  -. 1 739247E+0*CN+  . 87405830E-5*CN**2  -. 9883  I SCE- 1 C 

93* I *CN**3  _ _ 

94*  CTHR2=  -. II39E-4  -. I 296 1 I 0E+0*CN+  ,65869443E-5*CN**2  -. 7605028E- 10 

95*  l*CN**3 

96*  FDOTC2=“9288E‘-4  -.  I 40007 3E+0*Cf!+  _.696'444'|  6E-5*CH**2  -.  7556360E-1  0 

97*  __I*CH**3  

98*  XXX=( AMACH-O. 8 ) /0.2 

99*  GO  TO  72!  __  

100*  739  I F < AMACH- 1.1)740,740,741 
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SUST4  CONTINUED 


1.01  * . 740  CTHR I = 1 I 39E-4_r.J  296  1 1 QE+Q*CN+  .65869443E-5*CN**2  -.  76Q5C23E-  i 0 

1 02  * I *CN**3 

: 103.* FDOTC  i*__.92fi3E-4_=.,i.400073E+0*CN+  J>964  441  6H-5*CIJ**2-=.-7.556360E-.l 0- 

104*  1 *CN**3 

105*  CTHfi2=  -»11.93E-3--.J3I4.196E+0*CN+  .67603638£-5*CN**2  -.793421  5E-10 

106*  I *CN**3 

107* FDOTC2=- . J 7.68Er3_~.  962558 1 E- 1 *CN+  ..45264029E-5*CN**2  -.4 1 51 6 1 6E-I 0 

108*  1 *CN**3 

109* XXX=(AMACH-1.0)/0.1. 

.110*  721  CTHR=TERPI CXXX,CTHR2,CTHRI ) 

111* .FDOTC=TERPI  ( XXX.,  FD0TC2,  FDOTC  1 ) 

.112*  CO  TO  742 

113*  741  CTHR  = -.1I98E-3  -. 1 31 4I96E+0*CN+  .676C3638E-5*CN**2  79342 1 5E-1 0 

114*  l*CN**3 

1.15*  _.  ..  FDOTC  =-.  1 768E-3  -..962558  I E- 1 *CN+  ,45264029E-5*CN**2  -.4  1 51  6 1 6E-1  0 — 

.116*  1 *CN**3 

117*  .742  THR=T.W*DELTA*CTHR  _ 

.118*  FDC)TA=TM*FDOTC*DELTA*SQRTC THETA) 

119*  __  IF(THR.LT.0.0)THR=+ 1 .0  

120*  I FC ITD.GT. I 0 )THR=0. 0 

1 21  * I F ( ITD.GT.  1 0 ) FDOTA=0. 00  1 

1 22*  FISP=3600.*THR/FD0TA 

123*  IFCFISP.LT. 200. 00)FISP=200. 00  . 

124*  TITT2=3.06+0.00594*<ABS(RSPEED-75.0n**l  .71646 

125* TIT=TITT2*TT 

126*  WDOTAC=l 0.02-0, 0065* ( I 06.0-RSPEED)**l .87272 

1 27  * WDOT=WDOTAC*DELTA*TM/-SQRT  ( THETA ) 

I 28*  A0=0. 09  I 32*NDOT/C  R()*AMACH*SGRTCT ) 1/144.0 

129*  . . CAR=AO/AC/TM  — 

130*  CDA=0.0 

131*  IFCAMACH.GT. I .0)CDA=C 1 .51 7*(L0GC AMACH > ) **0. 41  * ( 1 .O-CAR) *AC/AREF 

132*  IFCCDA.LT. O.OOOOO: )CDA=0.0 

133* - CD0=CD()+CDA 

134*  THR()T=ASPEED 

135*  .THRS-THR  - 

136*  F ISPS=F ISP 

..  137*  799  CDAOFF=CPT+P)*AC*66.0/SQRTCTT)/( 1 .0+0. 2*AMACH**2. 0) **3. 0*SQRT< 5 .0* 

138*  1 ( (0 . 5*( PT+P ) /P) **( . 286 )- I .0  ) ) *( V-A*SQkT( 5.0*( CO. 5*( PT+P )/P ) ** ( .286 

139* 2)-l  .0) ) )/32. 1 74+TM/G/AREF*.. 

140*  IFCTHRS.GT.0.01 )CDAOFF=0.0 

. 1 41 *C  END  OF  SUSTAINER  DECK.  

142*  RETURN 

_ J 43* END  . 
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Appendix  I 
TRAJECTORY  EXAMPLES 


Five  example  trajectories  are  illustrated  in  this  appendix.  Each 
is  represented  by  a plot  and  demand  mode  output  format.  Each  plot  is 
annotated  so  that  the  user  can  correlate  the  program's  response  to  the 
input  data. 

The  output  column  labled  CDX  can  contain  either  the  total  zero- 
lift  drag  coefficient  (CDO) , or  the  inlet  spillage  drag  coefficient. 
Since  both  coefficients  will  be  based  on  AREF,  CDO  will  be  a number 
in  the  range  0.15  - 1.5,  and  spillage  drag  coefficient  will  be  con- 
siderably smaller  (0  - 0.1).  This  difference  in  magnitude  should 
provide  for  discrimination. 
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1234560000  0 0 15  0 I 5 

1200.  900.  1.58  15. ”.47  1.00'  5.0  5.00  . 
30.0  .90  10.0-10.040.00  1.40  100.  10.0  I. 
'STANDARD  ATMOSPHERE  ' 


0 0 
00  .0 
5010.00  I .60 


"TIME 

"ALT 

MACH 

RANGE" 

WEIGHT 

THRUST 

DRAG 

ALPHA 

L/D 

r PA  — 

CDX 

M 

SEC. 

FT. 

NO 

MM 

LB 

LB 

LB 

DEG 

DEG 

4.0 

5000. 

.13 

.0 

1 198.4 

1264. 1 

7. 

~ .0 

" ~.o 

~ "7  0 

723 

1 

8.0 

5000. 

.25 

.2 

1 196.7 

1 222.6 

26. 

.0 

.0 

.0 

.22 

1 

12.0 

"5000. 

.36 

' .4 

I 194.9 

“T  21 I .0 

""  56. 

.0 

.0" 

— .0  - 

7 22 

1 

TAKE-OFF  ROLL 

REQUIRED  3797.  FEET. 

17.3 

4988. 

.50 

.8 

1 192.5 

1 230.9 

532." 

15.0 

4.0 

1 .3" 

722 

2 

27.3 

7003. 

.60 

1 .7 

1187.7 

1215.6 

214. 

5.2 

4.7 

29.0 

.22 

2 

— 37.3 

10529. 

.70 

2.7 

1183.0 

"1161.0 

226. 

4. 1 

4.2 

30.0 

.22 

2 

47.3 

1 4509 . 

.79 

3.9 

1 1 78.7 

1096. 1 

240. 

3.8 

4.0 

30.0 

.23 

2 

57.3 

1 8841 . 

“86 

5.T 

(174.7  1017.  r 

“244." 

” 3.8 

"3.9 

"30.0" 

.23 

2 

67.3 

2341 7. 

.91 

6.4 

1171.0 

922.7 

237. 

4.  1 

4.0 

30.0 

.24 

2 

77.3 

28117. 

.‘/'4 

' '7.3 

1167.8 

' 814.2 

235. 

4.7 

4.0 

30.0 

.25 

' 2 

87.3 

32677. 

.95 

9.1 

1 165.0 

706.3 

230. 

5.6 

4.1 

28.5 

.27 

2 

~ '97.3 

36966. 

.93 

10.5 

1162.6 

600.6 

208. 

7.1 

4.5 

27.2 

.26 

2 

107.3 

39926. 

.95 

1 1 .9 

1 160.6 

51  7.3 

109. 

.0 

.0 

9.6 

.28 

2 

111.3 

40292. 

.98 

12.5 

1159.9 

" "523.4 

1 47. 

.0 

.0 

1 .7 

.35 

3 

1 15.3 

40228. 

1 .01 

13.1 

1159.1 

526.4 

368. 

8.7 

3.1 

-1  .8 

.00 

3 

1 27.3 

40036 . 

1 .06 

15.1 

1156.8 

" 547.2 

407. 

7.8 

2.7 

-.3 

.00 

3 

147.3 

40002. 

1.15 

18.6 

1152.6 

585.5 

441  . 

6.5 

2.5 

-.0 

.01 

3 

l”67.3 

40000. 

1.23 

22.4" 

1148.1 

622.1 

468. 

5.6 

2.3 

-.0 

.01 

3 

187.3 

40000. 

1.32 

26.5 

1143.4 

660.0 

500. 

4.9 

2.2 

-.0 

.01 

3 

207.3 

40000. 

1.40 

30.3 

11 38.4 

627.5" 

537. 

4.3 

2.0 

' -.0 

.01 

'3 

227.3 

40000 . 

1 . 41 

35.3 

1 134.3 

519.5 

543. 

4.3 

2.0 

-.0 

.01 

3 

247.3 

40000. 

1 . 40 

“39.8 

1 130.6 

521.9 

536. 

4.3 

2.0 

.0 

.01 

3 

267.3 

40000. 

1.40 

44.2 

1126.8 

540.2 

534. 

4.3 

2.0 

.0 

.01 

3 

287.3 

40000. 

1.40 

48.7 

1122.9 

538.8 

" 535. 

4.3 

2.0 

-.0 

.01 

3 

307.3 

40000. 

1.40 

53.2 

1 119.0 

534.6 

534. 

4.3 

2.0 

-.0 

.01 

3 

327.3 

40000. 

1 .40 

57.6 

11 15.2 

533.9“ 

533. 

4.3 

2.0 

-.0 

.01 

3 

347.3 

40000. 

1 . 40 

62.1 

1 1.1  1 .3 

533.7 

532. 

4.3 

2.0 

-.0 

.01 

3 

367.3 

40000. 

1 . 40 

66.5" 

I 107.5 

532.9 

"531  . 

4.2 

2.0 

-.0 

.01 

3 

387.3 

40000. 

1.40 

71  .0 

1103.6 

532.0 

531  . 

4.2 

2.0 

-.0 

.01 

3 

“ 407.3 

40000. 

1 .40 

75.5 

1099.8 

"531.2 

530. 

4.2 

2.0 

-.0 

.01 

3 

427.3 

40000. 

1 .40 

79.9 

1096.0 

530.4 

529. 

4.2 

2.0 

-.0 

.01 

3 

447.3 

"40000. 

1 . 40 

*84.4" 

1092. 1 

"529.7 

528. 

4.2 

2.0 

“-.0 

.01 

3 

467.3 

40000. 

1 .40 

88.9 

1 088. 3 

528.9 

527. 

4.2 

2.0 

-.0 

.01 

3 

487.3 

40000. 

1 .'40 

93.3 

1084.5 

528.  1 

527. 

4.2 

2.0 

-.0 

.01 

3 

507.3 

40000. 

1.40 

97.3 

1080.7 

527.3 

526. 

4.  I 

2.0 

-.0 

.01 

3 

“519.3 

39939. 

1 .41 

100.5 

1078.4 

699.3 

424. 

.0 

.0 

-2.6 

.01 

4 

523.3 

39424. 

1 . 46 

101  .4 

1077.3 

732.5 

457. 

.0 

.0 

-7.9 

.01 

4 

“527.3 

38476. 

1 .51 

102.3 

1076.2" 

710.0 

"589  .“ 

"""  3.2 

"1.7 

-10.0 

.or 

" “4 

547.3 

33280. 

1.53 

107.1 

1072.3 

459.0 

708. 

2.6 

1.5 

-10.0 

.01 

4 

~ 567.3 

28125. 

1.49 

1 12.0 

1068. 1 

“ 675.1 

814. 

~ 2. 1 

1 .3 

-10.0 

.01 

4 

587.3 

22898. 

1.49 

1 16.8 

1062.2 

859.1 

984. 

1 .7 

1 .0 

-10.0 

.01 

4 

607.3 

17547. 

1 . 49 

121  .8 

1054.7 

1063.2 

1193. 

1.3 

' .8 

-10. 0 

.01 

4 

627.3 

1 2081 . 

1 . 49 

1 26.9 

1 045.1 

1316.2 

1 444. 

1.1 

.7 

-10.0 

.01 

4 

“635.2 

10018. 

1.48 

129.0“ 

1040.4 — 1782.0 

2162. 

“5.0 

“2.3 

“ -3.9— 

.01 

"5 

637.2 

9943. 

1.48 

129.5 

1033.3 

1 810.8 

1605. 

1.5 

.9 

.0 

.01 

5 

“640.4 

9966. 

1.50 

130.3 

1036.3 

"1826.5 

1603. 

.9 

“ .6 

.3 

.01 

5 

650.4 

9996. 

1.55 

133.0 

1028. 4 

1857.7 

1704. 

.9 

.6 

.0 

.01 

5 

"660.4 

10000. 

1.59 

135.8 

1020.4 

1860.3 

1770. 

.8 

.6 

“ .0 

.01 

5 

670.4 

1 0000. 

1.60 

138.7 

1012.6 

1 791 .2 

1790. 

.8 

.6 

.0 

.01 

5 

680.4 

1 0000. 

1.60 

141 .5  1004.9 

~ 1 786.2 

7 7 88 

. 8 

75 

. 0 

.01 

~ 5' 
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NWC  TM  3540 


— TIME 

ALT 

MACH 

RANGE 

WEIGHT 

THRUST  DRAG  ALPHA  L/D 

' FPA 

_ CDX  ' 

M 

690.4 

1 OOOO. 

1 .60 

1 44.3 

997.3 

1787.6  1737. 

.8  .5 

.0 

.01 

5 

— 700.4 

1 0000 . 

1 . 60 

147.2 

989.7 

~ 1 737. 3 1787. 

—.8  1 5 

• 0 

•Of 

5 

710.4 

10000. 

1.60 

150.0 

982.0 

1786.9  1787. 

.8  .5 

.0 

.01 

5 

—720.4 

1 0000 . 

1 .60 

152.9 

974.4 

1786.6  1786. 

.8  -.5 

— .0 

—.01 

5 

730.4 

10000. 

1 .60 

1 55.7 

966.8 

1786.3  1786. 

.8  .5 

.0 

.01 

5 

— 740.4 

1 0000. 

1.60 

158.5 

959.  1 

1786.0  1786. 

.8'— .5  ‘ 

- .0 

—.01  - 

5 

750.4 

1 0000 . 

1.60 

161.4 

951.5 

1735.7  1786. 

.8  .5 

.0 

.01 

5 

7 60 . 4 ' 

1 0000 . 

1 .60 

164.2 

943.9 

' 1 785.3  1 785. 

.8  .5  “ 

.0  .XM 

5 

770.4 

1 0000. 

1.60 

167.0 

936.3 

1785.0  1785. 

.7  .5 

.0 

.01 

5 

—700.4 

10000. 

1 . 60 

169.9 

928.5 

1 784.7  1 735. 

. 7 .5 

— .0' 

• 0 1 

5 

790.4 

1 0000. 

1 .60 

1 72.7 

921.0 

1734.4  1734. 

.7  .5 

.0 

.01 

5 

— 800.4 

1 0000 . 

1.60 

175.5 

913.4 

1784.1  1784. 

.7  -.5- 

— -o 

—.01  - 

5 

810.4 

10000. 

1.60 

1 78.4 

905.8 

1783.8  1784. 

.7  .5 

.0 

.01 

5 

— 820.4 

~ 9944. 

I .47 

] 81  .2 

— 899.9 

. 0‘ 1 752 l 

• .o  ;o— 

-Z7T 

.01 

6 

830.4 

7999. 

1.08 

183.4 

899.9 

.0  972. 

.0  .0  - 

•15.9 

.01 

6 

—840.4 

3701  . 

.96 

185.0 

899.9 

.0  513. 

.0  .0  - 

■31.6 

—.01 

6 

“S"  S — 

— S 

SSSS 

S"~ 

55555  55555  55555 

S S~  ! 

s s 
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s 

S 

SSSS 
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"SSS 

sssss 
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5- 

S "S  ~ S SSSS 

" S ~ ! 

s s 

s s 

S S 

s 

S 

S S s s 

S ! 

-S  S - 

s s 

SSSS 

sssss  sssss  sssss  sssss 

S 0 

— 846.4 

0. 

.94 

185.9 

899.9 

.0  539. 

.0  .0  - 

40.6 

-01- 

6 

~IT  HAS  BEEN  A PLEASURE  SERVING  YOU.  I HOPE  YOU  ENJOY  YOUR  "T ‘FRAMES  OF'FR-80 
UT.  


1-4 


77 


EXAMPLES 


RANGE 


NWC  TM  3540 


1638540000  0 0 15  0 I 5 0 0 

~l  200.  900.  1.58  15.  .47  '1.00  5.0  5.00  .00  “.0  '"  ‘ 

20.0  .90  10.0-10.040.00  1.40  JOO.-IO.O  1.5010.00  1.60 

-STANDARD  ATMOSPHERE  


-TIME 

- ALT  - 

MACH 

RANGE 

WEIGHT 

THRUST 

DRAG 

ALPHA 

L/D 

- FPA 

CDX 

SEC. 

FT. 

NO 

NM 

LB 

LB 

LB 

DEG 

DEG 

-4.0 

5000. 

-.13 

.0 

1198.4 

1264. 1 

7. 

" .0 

.0 

.0 

.23 

"1 

8.0 

5000. 

.25 

.2 

1196.7 

1 222.6 

26. 

.0 

.0 

.0 

.22 

1 

12.0 

5000. 

. 36 

“ .4 

1 194.9 

1211.0 

56 . 

-.0 

.0 

.0 

.22 

"I 

TAKE-OFF  ROLL 

REQUIRED  3797.  FEET. 

— 17.3 

5059. 

.49 

.8 

1192.5 

1227.2 

504. 

1 5.0 

4.0 

" 6.3 

.22 

6 _ 

27.3 

7449. 

.55 

1 .7 

1187.7 

1247.3 

319. 

8.9 

4.9 

39.9 

.22 

6 

37.3 

.1  1862. 

.59 

2.4 

1183.1 

1187.4 

190. 

4.5 

4.3 

47.6 

.23 

6 

47.3 

1 6465. 

.62 

3.2 

1 179.2 

1067.2 

1 76. 

4.3 

4.1 

43.4 

.23 

6 

57.3 

20742. 

.64 

- "4.0 

-1175.7 

946.6 

“173. 

— 5.1 

4.4 

38.0 

.24 

6 

67.3 

24536. 

.65 

4.9 

1172.7 

840.5 

174. 

6.  1 

4.6 

33.0 

.25 

6 

—77.3 

27979. 

.66 

5.9 

3170.0 

750.6 

1 78. 

7.3 

4.7 

'28.5 

.25 

6 

87.3 

30935. 

.67 

6.8 

1 167.6 

675.7 

186. 

8.6 

4.7 

24.5 

.26 

6 

97.3 

33481 . 

.67 

7.9 

1165.4 

'61 4.0 

196. 

1 0. 1 

4.5 

21  .0 

.26 

6 ' " 

107.3 

35651 . 

.67 

8.9 

1 163.4 

563.5 

210. 

1 1 .6 

4.3 

17.9 

.26 

6 

— 117.3 

37479. 

‘.67 

-"9.9" 

"1 161  .6 

~ 522.5' 

224. 

" 1 3"  1 

4.1 

15.0 

.26 

6 

127.3 

38978. 

-56 

10.9 

1159.9 

478.0 

239. 

14.7 

3.9 

12.2 

.27 

6 

— 1 36 .6 

40046. 

-.66 

11  .9 

1158.5 

444.5 

53. 

.0 

.0 

7.6 

.27 

"3 

140.6 

40146. 

.69 

12.3 

1157.9 

430.0 

239. 

1 5.0 

3.3 

-1  .3 

.27 

3 

— 144.6 

40061 . 

.72 

12.8 

1157.3 

433.3 

255. 

1 5.0 

3.8 

-2.2 

.27 

3 

158.6 

39701 . 

.79 

14.4 

1155.1 

456.6 

311. 

15.0 

3.8 

-.7 

.26 

3 

1"78.6 

39957. 

.91 

~f7.r 

-JI5I  .7 

“ 502.6 

229. 

9.8 

4.5 

.4 

.26 

—3 

193.6 

3 9998. 

1 .03 

20.2 

1147.9 

541.6 

366. 

7.9 

2.9 

.0 

.45 

3 

— 218.6 

40000 . 

1.12 

23.7 

1 1 43.8 

‘ 576.2 

425. 

‘ 6.7 

2.5 

.0 

.48 

_3 

238.6 

40000. 

1 .21 

27.4 

1 139.5 

612.9 

454. 

5.8 

2.4 

.0 

.48 

3 

— 258.6 

40000. 

1.30 

-31  .4 

1134.7 

652.6 

“486. 

“5.0 

2.2 

— .0 

.48 

”3  — 

278.6 

4C000 . 

1.40 

35.7 

1 129.8 

64  1 .9 

524. 

4.3 

2. 1 

-.0 

.48 

3 

—298.6 

40000. 

1.42 

40.2 

1 125.6 

“ 516.5 

536. 

4.2 

~ 2.0 

-- .0 

”.48 

"73  — 

318.6 

40000 . 

1 .40 

44.7 

1 122.0 

509.0 

523. 

4.3 

2.1 

.0 

.48 

3 

- '338.6 

40000. 

I . 40 

49. 1 

1118.2 

‘530.6 

525. 

4.3 

2.1 

-.0 

".48 

3 

358.6 

40000. 

1 . 40 

53.6 

1114.4 

530.8 

526. 

4.3 

2.0 

-.0 

.48 

3 

— 378.6 

40000. 

1 .40 

58.0 

1110.6 

" 526.1 

525. 

4.3 

2.0 

-.0 

.48 

3 — 

398.6 

400C0. 

1 .40 

62.5 

1106.8 

525.2 

524. 

4.2 

2.0 

-.0 

.48 

3 

—418.6 

40000 . 

1.40 

67.0 

“1  103.0 

—525.1 

524. 

4.2 

2.0 

-.0 

.48 

3 

438.6 

40000. 

1 .40 

71  .4 

1099.3 

524.3 

523. 

4.2 

2.0 

-.0 

.43 

3 

'458.6 

40000. 

1 . 4Q 

75.9 

1095.5 

- 523.4 

522. 

4.2 

2.0 

-.0 

.48 

3 

478.6 

40000. 

1 .40 

30.3 

1091.7 

522.6 

521  . 

4.2 

2.0 

-.0 

.48 

3 

'498 . 6 

40000. 

1.40 

84.8 

1088.0 

521.9 

521  . 

" 4.2 

2.0 

-.0 

.48 

”3 

518.6 

40000 . 

1 .40 

89.3 

1084.2 

521.1 

520. 

4.2 

2.0 

-.0 

.48 

3 

—538.6 

40000. 

1 . 40 

-"93.7-1080.5 

— 520.4 

" 519. 

— 4.1' 

— 2.0 

— _vcr 

—.48 

3 

558.6 

40000. 

1 .<0 

98.2 

1076.7 

51  9.6 

51  8. 

4.1 

2.0 

-.0 

.48 

3 

“ 572.8 

39487. 

1.45 

101  .4 

1073.8 

670.7 

437. 

.0 

.0 

-7.5 

.48 

8 

582.8 

36558. 

1.54 

103.7 

1071.6 

487.0 

629. 

2.9 

1 .6 

-12.9 

.47 

8 

"592.8 

33272. 

1.51 

106.1 

1070. 1 

- 393.1 

693. 

— 2.6 

1 .5 

-12.6 

.46 

8 

602.8 

30087. 

1 . 49 

103.5 

1068.4 

568.8 

754. 

2.3 

1 .4 

-12.3 

.46 

8 

— 612.8 

26956. 

1 .49 

‘110.9 

1065.9 

— 722. 1 

— 845 . 

2.0 

‘‘  1 .2 

-12.0 

.46 

"8 

622.8 

23847. 

1.50 

1 13.3 

1063.0 

805.6 

950. 

1.7 

1 . 1 

-11  .7 

.45 

8 

' 632.8 

20758. 

1.50 

115.8 

1059.7 

907.6 

1060. 

1 .5 

1 .0 

-11  .5 

.45 

3 — 

642.8 

1 7686. 

1 .49 

1 18.3 

1055.9 

1036.7 

1182. 

1 .3 

.9 

-II  .3 

.44 

8 

' 652.8 

1 4625. 

1.50 

120.8 

1051.5 

1175.2 

1317. 

1.2 

.8 

-11.2 

.44 

8 

662.8 

11567. 

1.50 

123.4 

1046.4 

1324.7 

1464  . 

1 .0 

.7 

-1  1 .0 

.43 

8 

-'  "668.0 

1 0067. 

1.49 

124.7 

1043.3 

‘ 1797.4 

2151  . 

4.9 

—2.4 

-5.5 

_.43 
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NWC  TM  3540 


~ TIME 

ALT 

MACH 

RANGE 

HEIGHT 

THRUST 

DRAG  ALPHA 

L/D 

FPA 

CDX 

M 

670.0 

9947. 

1 .48 

125.3 

1041.7 

1811.4 

1606. 

1 .8 

1 .2 

-.3 

.43 

5 

— 672.6 

"’9960. 

1.50' 

125.9 

1039. T 

"'1827. 1 

I 579 . 

.9' 

. 6 

' /3" 

.43 

5 

682.6 

9996. 

1 .56 

128.7 

1031.7 

1862.6 

1685. 

.9 

.6 

.0 

.42 

5 

~ 692.6 

10000. 

1 .60 

131  .5 

1023.8 

1791.6 

1747. 

.8 

.6 

.0 

~742 

5 

702.6 

1 0000 . 

1.60 

134.3 

1016.3 

1 745.8 

1 751  . 

.8 

.6 

.0 

.42 

5 

"■‘712.6 

1 0000. 

1 

137. 1 

1003.9 

1748.6 

1749. 

.8 

.6 

-.0 

~ .42 

5 ' 

722.6 

1 0000 . 

I .60 

1 40.0 

1001.5 

1749.1 

1 749. 

.8 

.6 

.0 

.42 

5 

732.6 

10000. 

1.60 

142.8 

994.0 

"1748.7 

1749. 

.8 

“'.6" 

—.0 

.42 

5 ' 

742.6 

1 0000  . 

I .60 

145.7 

986.6 

1 748.4 

17  48. 

.8 

.6 

.0 

.42 

5 

752.6 

1 0000 . 

1 .60 

1 48.5 

979.  1 

1748. 1 

1748. 

.8 

.5 

.0 

■'.42 
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.24" 

3 

328.8 

20000. 

.96 

40.9 

1137.9 

416.2 

394.4 

3.0 

2.8 

-.0 

.26 

3 

“348.8 

20000. 

-'.96 

44.2 

1 134.5 

"417.4 

415.2 

“ 3.0 

2.7" 

““  -.0 

.28" 

3" 

368.8 

20000. 

.96 

47.5 

1131.1 

41  7.6 

416.7 

3.0 

2.7 

-.0 

.28 

3 

— 383.3 

20000 .' 

'“.96" 

"50 . 8 ' 

1127.-7 

417.6" 

416.9 

"2.9 

2.7' 

-.0 

“28“ 

3“ 

403.8 

20000. 

.96 

54.1 

1 124.3 

41  7.6 

416.9 

2.9 

2.6 

-.0 

.28 

3 

— 428.8 

20000. 

.96 

57.4 

1120.9 

— 417.6 

416.9 

2.9 

2.6 

“ -.0 

.28  "3  “ 

448.8 

20000. 

.96 

60.7 

1 1 1 7.5 

41  7.7 

417.0 

2.9 

2.6 

-.0 

.28 

3 

— 468.8 

20000. 

.97 

"64.0- 

iiu.r 

417.7“  - 

417.0 

2.9 

“2.5 

-70 

“ .28  “ 

~3 

488.8 

20000. 

.97 

67.3 

1110.7 

417.7 

417.0 

2.9 

2.6 

-.0 

.28 

3 

— 508.8 

'20000. 

.97  7075 

71 07.3' 

417.7 

417.0 

279  ' 

2.6 

-.0 

“ ."23“ 

5 

528.8 

20000. 

.97 

73.8 

1 103.9 

41  7.7 

417.1 

2.9 

2.6 

-.0 

.28 

3 

— 548.8 

20000. 

.97 

77.1 

1 100.5 

' “ 417.8 

417.1 

“2.9 

2.6 

-.0 

728 

3 “ 

568.8 

20000. 

.97 

30.4 

1097.  1 

41  7.8 

417.1 

2.9 

2.6 

-.0 

.28 

3 

“588.8 

20000. 

.97 

33.7 

1093.7 

417.8 

417.1 

2.8 

2.6 

--.0 

" .28  " 

3'  “ 

608.8 

20000. 

.97 

37.0 

1090.3 

417.8 

417.2 

2.8 

2.6 

-.0 

.28 

3 

628.8 

20000.' 

.97- 

90.3  1 086.9" 

"41^79 

“417.2' 

“ 278 

“2.6 

-.0' 

3 
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TIME 

ALT 

MACH 

RANGE 

WEIGHT 

THRUST 

DRAG 

ALPHA 

L/D 

FPA 

CDX 

M 

—648.8 

20000. 

.97 

-93.6 

(083.5 

41 7.9 

417.2 

2.8  - 

2.5 

---.0 

728 

—3 

668.8 

20000. 

.97 

96.9 

1 080. 1 

417.9 

417.3 

2.8 

2.5 

-.0 

.28 

3 

— 688.8 

20000. 

.97 

100.2 

1076.7 

— 417.9 

417.3 

2.8  - 

2.5 

- “-.0 

—28 

3 

708.8 

20000. 

.97 

103.5 

1073.3 

41  7.9 

417.3 

2.8 

2.5 

-.0 

.28 

3 

728.8 

20000.- 

.97 

1 06 . 8 

1 069.9- 

— 41 8.0 

-417.3 

2.8 

2;5 

— -.0 

— T28 

-3 

748.8 

20000. 

.97 

110.1 

1066.5 

41  3.0 

417.4 

2.8 

2.5 

-.0 

.28 

3 

— 768.8 

20000. 

.97 

1 13.4 

1063.1 

— 41 8.0 

417.4 

2.8 

2.5 

“ -.0—728 

_ -3  

788.8 

20000. 

.97 

116.7 

1059.7 

413.0 

417.4 

2.8 

2.5 

-.0 

.28 

3 

- 808.8 

20000. 

.97 

120.0 

1056.3 

418.1 

417.4 

2.7 

2.5 

.0 

-.28 

3 

828.8 

20000. 

.97 

123.3 

1052.9 

418.1 

417.5 

2.7 

2.5 

-.0 

.28 

3 

— 848.8 

20000. 

. 97 

" I 26  • 6 

1 049.5 

418.1 

417.5 

2.7  - 

2.5 

-.IT 

3 

868.8 

2C000 . 

.97 

1 29.9 

1 046. 1 

41  8.  1 

417.5 

2.7 

2.5 

-.0 

.28 

3 

— 888.8 

20000. 

.97 

133.2 

1042.6 

—418.1 

417.5 

2.7 

2.4 

-.0 

~ .28 

3 

908.8 

20000. 

.97 

1 36.5 

1039.2 

41  8.2 

4 17.6 

2.7 

2.4 

-.0 

.28 

3 

—928.8 

20000. 

.97 

139.8 

1035.8 

~ -418.2 

417.6 

2.7 

2.4 

- -.0 

~29 

-3 

948.8 

20000. 

.97 

1 43.1 

1032.4 

418.2 

417.6 

2.7 

2.4 

-.0 

.29 

3 

— 968.8 

20000. 

.97 

I 46.4- 

1029.0 

— 418.2 

417.6 

" 2.7 

2.4 

-.0 

729" 

-3 

988.8 

20000. 

.97 

1 49.7 

1025.6 

418.2 

417.7 

2.7 

2.4 

-.0 

.29 

3 

-1008.8 

1 8660. 

- .90 

152.8 

1024.3 

— 104.9 

-305.5- 

3.P 

3.3 

-5.0 

"-.23 

"9 

1028.8 

1 7100. 

.83 

155.8 

1023.4 

107.5 

286.9 

3.4 

3.5 

-5.0 

.23 

9 

“1048.8 

1 5645. 

.77 

158.5 

1022.5 

- 1 10.3 

273.8 

3.7 

3.7 

* -5.0 

“ .23 

9 - 

1068.8 

1 4278. 

. 72 

161 . 1 

1021.6 

113.3 

263.8 

4 . C 

3.3 

-5.0 

.23 

9 

-1088.8 

1 2932. 

.63 

"1  63. 5~1  020.6" 

115.6 

~ 256.0 

4.2  ‘ 

3.9 

-5.0' 

““23 

9 ' 

1108.8 

.1  1747. 

.65 

165.9 

1019.6 

1 1 7.2 

249.9 

4.4 

4.0 

-5.0 

.23 

9 

— 1128.8 

1 0564 

' .62 

168. 1 

101 3.5 

T18.3 

“244.9 

4.6  ‘ 

4.1 

“ -5.0 

“723 

-9 

1 1 48.8 

9425. 

.60 

1 70.2 

1017.4 

1 1 9.4 

240.8 

4.8 

4.2 

-5.0 

.22 

9 

“ TT68.8 

83251 

-.53 

172.3 

101 6.3 

— 122.8 

237 .6 

“4.9  - 

4.2 

'-5.0 

“7  22 

9 

1 1 88.3 

7257. 

. 56 

1 74.3 

l 01 5. 1 

125.4 

235.2 

5.0 

4.3 

-5.0 

.22 

9 

— T208.8 

' ’6215. 

'.54 

176.3 

1 01 3.9 

— 127.4 

“233/3 

“5.0 

4.3' 

“ -5.0' 

“122 

9 

1228.3 

5195. 

.53 

1 78.2 

1012.6 

129.0 

231  .8 

5.  1 

4.3 

-5.0 

.22 

9 

1 248.8 

4194. 

.52 

1 80.1 

1011.3 

— 130.1 

230.5 

5.  1 

4.3 

-5.0 

.22 

9 

1268.8 

3208. 

.51 

1 81  .9 

1009.9 

131  .0 

229.5 

5. 1 

4.3 

-5.0 

.22 

9 

- 1288.8 

2236. 

.50 

183.7 

1003.5 

131.5 

228.5 

5.0 

4.3 

-5.0 

.22 

9 

1307.6 

1 335. 

.49 

1 85.4 

1007. 1 

131.8 

246.3 

5.5 

4.4 

-4.9 

.22 

9 

-1311 .6 

1 1 74. 

.49 

1 85.8 

1006.8 

131.7 

246.5 

5.5 

4.4 

-3.5 

.22 

9 

1315.6 

1068. 

.48 

186.1 

1006.5 

133.8 

2 45.3 

5.6 

4.5 

-2.1 

.22 

9 

-1319.6 

1015. 

.47 

186.5 

1006.2 

136.7 

244.0 

5.8 

4.5 

-.7 

.22 

9 

1323.6 

1007. 

. 46 

1 86.8 

1 005.9 

139.9 

223.5 

5.5 

4.4 

-.1 

.22 

9 

-1327.9 

1006. 

. 45 

1 87.2 

1005.6 

142.9 

22!  .9 

5.8 

4.5 

-.0 

.22 

9 

TOUCHDOWN  AT  FLIGHT  PATH  ANGLE  OF  -.04  DEGREES 


1345.3  1000 

. — .27  1 88.3  1005.0 — 

.0“ 

- ■4 1.9  “7  0 “70  .0  —.731“ 

9 

LANDING  COMPLETED  TO 
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STOP  IN 

7159. 
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7831000000  1 0 15 

— 1 763.  900.  1.58  15.  .47 

20.0  .90  30.0-10.015.00 

STANDARD  ATMOSPHERE 

6 

- 1 .00- 
1.40 

1 5 
-5.0  2.0C 
80.  10. C 

0 0 
_.00  40.0 

1 .3010.00 

.60 



— 

_ TIME 

ALT  MACH 

RANGE 

WEIGHT 

__  THRUST 

DRAG 

ALPHA  L/D 

FPA 

CDX 

M 

SEC. 

FT.  NO 

NM 

LB 

LB 

LB 

DEG 

DEG 

I .0 

-2043.  - .23 

.0 

1707.2 

13575.9 

80.5 

9.7 

—3.9 

-20.4 

.36 

a - 

2.0 

2182.  .43 

.1 

1658.6 

12105.7 

222.9 

7.  1 

3.8 

22.9 

.36 

8 

3.0 

2419.  .62 

.2 

1614.8 

.11208.8 

386.4 

5.4 

3.4 

24.6 

.35 

8 

4.0 

2757.  .80 

.3 

1574.5 

10591.9 

556.7 

4.  1 

2.9 

26.0 

.35 

8 

5.0 

3194.  .97 

.4 

1537.5 

9297.3 

833.6 

3.0 

2.0 

27.1 

.39 

8 

6.0 

3717.  1.07 

.6 

1513.7 

4964.2 

1425.0 

2.3 

1.1 

27.7 

.02 

S 

7.0 

4277.  1.10 

.8- 

1505.4 

2439. 1 

1490.6 

2.0 

.9 

- 28.0 

.02 

8-- 

1 1 .6 

6891.  1.10 

1 .6 

1396. 1 

161 8.3 

1056. 1 

1.8 

1 . 1 

28.2 

.02 

8 

21  .6 

.11922.  1.11 

3.3 

1389.8 

1434.5 

863.9 

.0 

.0 

19.3 

.02 

8 

31  .6 

14733.  1.16 

5.3 

1334.  1 

1364.4 

844.  I 

.0 

.0 

6.6 

.02 

3 

37.2 

15045.  1.22 

6.4 

1 380.9 

1396.8 

983.1 

2.2 

1 .3 

-.4 

.02 

3 

41  .2 

15019.  1.25 

7.3 

1378.6 

1425.3 

1043.0 

2.2 

1 .3 

-.2 

.02 

3 

48.8 

1 5004 . 1.31 

- 9.0 

1374.0 

1 475.4 

1131.6 

2.0 

1 .2 

-.0 

.03 

3 

68.8 

15000.  1.40 

13.7 

1362.2 

1327.8 

1284.3 

1.7 

1 .0 

-.0 

.03 

3 

08.8 

14656.  1.39 

18.6 

1352.1 

966.  1 

1629.1 

4.6 

2.  1 

-4.  1 

.03 

10 

108.8 

13679.  1.30 

23.3 

1344.0 

1072. 1 

1 1 70 . 7 

1.9 

1 .1 

-1  .8 

.03 

10 

— 128.8 

12827.  1.30 

27.8 

1334.9 

J 1 6 1 . 4 

1 1 88 . 8 

1.8 

1 . 1 

-1.8 

.03 

10 

148.8 

11992.  1.30 

32.4 

1325.3 

.1194.  1 

1 223.6 

1.7 

1 . 1 

1 .7 

.03 

10 

168.8 

1 1 183.  1 .30 

-37.0-1315.3 

-1226.8 

1256.0 

1.7 

1 .0 

— 1.6 

.03 

10 

1 88.8 

10400.  1.30 

41  .5 

1 305.1 

1260.1 

1287.9 

1 .6 

1 .0 

-1  .6 

.03 

10 

— 208.8 

-9644.  1.30 

-46. 1 

1294.5 

-1293.1 

1319.6 

1 .6 

-1.0 

1 .5 

.03 

10 

228.8 

8914.  1.30 

50.8 

1 283.7 

1325.6 

1 350.9 

1.5 

.9 

-1  .5 

.03 

10 

—248.8 

8212.  i-.  30 

-55.4 

1272.5 

—1357.6 

1381.7 

1.4 

.9 

- -1  .4 

.03 

10 

268.8 

7537.  1.30 

60.0 

1260.9 

1389.1 

1411.9 

1 .4 

.9 

-1  .3 

.03 

10 

— 288.8 

6890.  1.30 

64.7- 

1 249.  1 

—141  9.9 

1441.4 

1 .4 

— .8 

— 1 .3 

.03 

10 

308.8 

6271.  1.30 

69.3 

1237.0 

1 449.9 

1470.3 

1 .3 

.8 

-1.2 

.03 

10 

—328.8 

5680.  1.30 

-74.0 

1 224.6 

1479.1 

1 498.2 

1.3 

.8 

--1  .2 

.03 

10 

348.8 

5118.  1.30 

78.7 

1211.9 

1507.4 

1525.3 

1.2 

.8 

-1  . 1 

.03 

10 

- 368.8 

4584.  1.30 

83.4- 

1198.9 

1534.7 

1551  .4 

1 .2 

.8 

-1  .0 

.03 

10 

388.8 

4080.  1.30 

88. 1 

1185.6 

1 560.9 

1576.5 

1.2 

.7 

-1  .0 

.03 

10 

— 408.8 

- 3605.-1 .30 

-92.8- 

1172.1- 

—1586.0 

1 600 . 4 

1 . 1 

.7 

-.9 

.03 

10- 

428.8 

3159.  1.30 

97.5 

1158.4 

1609.8 

1623.2 

1 . 1 

.7 

-.9 

.03 

10 

— 448.8 

2742.  1.30 

102.2  - 

1144.4 

1632.5 

1644.7 

1 . 1 

- .7 

— -.8 

.03 

10 

468.8 

2355.  1.30 

107.0 

1130.2 

1653.7 

1664.9 

1 .0 

.7 

-.7 

.03 

10 

—488.8 

1997.  1.30 

111.7 

1115.7 

-1673.6 

1633.8 

1 .0 

.6 

— - .7 

.03 

10 

508.8 

1670.  1.30 

116.5 

1101.1 

1692. 1 

1701  .2 

1 .0 

.6 

-.6 

.03 

10 

— 528.8 

- 1371 .-1,30 

121  .2 

1086.3- 

-1-709.  1 

1717.2 

—1.0 

— .6- 

-.6 

-.03- 

10- 

548.8 

1102.  1.30 

1 26.0 

1071.3 

1724.6 

1731  .7 

.9 

.6 

-.5 

.03 

10 

— 568.8 

- 863.  1.30 

1 30.7 

1056.2 

-1738.5 

17  44.7- 

— .9 

.6 

-.4 

-.03 

10 

588.8 

654.  1.30 

135.5 

1041.0 

1750.9 

1756. 1 

.9 

.6 

-.4 

.03 

10 

— 608.8 

— 473. -1 .30 

1 40.2 

1025.6 

-1761 .6 

1765.9 

- .9 

--.6 

— -.3 

—.03 

10 

623.8 

322.  1.30 

1 45.0 

1010.  1 

1770.7 

1774.2 

.9 

.6 

-.3 

.03 

10 

—648.8 

201  . -L.30 

1-49  ,8  - 

-994.5- — 1-773.2 

1780.8 

—.8 

.5 

-.2 

— 03- 

-10- 

668.8 

108.  1.30 

1 54.6 

978.9 

1 784.0 

1785.8 

.8 

.5 

-.2 

.03 

10 

—674.4 

92. -1 .31 

1 55.9- 

-974.3 

-2090. 4 

1821  .7 

-.8 

- .5 

0 

.03 

5 

680.6 

98.  1.35 

157.4 

968.5 

2130.8 

1 940.0 

.7 

.5 

.0 

.03 

5 

- 690.6 

—100.  1.40 

1 59.9 

-958.9 

-2170.5 

2068.4 

7 

-.5- 

.0 

.03 

5 

700.6 

100.  1.42 

162.5 

949.  1 

2191 . 1 

2132.4 

.7 

.4 

.0 

.03 

5 

710.6 

. i no.  i-4i 

1 65.2 

939.2- 

-2202.7- 

2169.1- 

,6 

-.0 

—03 

5 

720^6 

100.  1.44 

167^8 

929^3 

2209.2 

2190.1 

.6 

'.4 

-'.6 

!o3 

5 
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TIME 

ALT 

MACH  RANGE 

weight  thrust  drag 

ALPHA  L/D 

FPA 

CDX 

M 

730.6 

1 00. 

1 . 45  1J0.4 

91 9.4  .221 2.8  2202.0 

.6  .4 

-.0 

.-  .03. 

5-  ...  - 

740.6 

100. 

1.45  173.1 

909.5  221 4.9  2208.7 

.6  .4 

-.0 

.03 

5 

750.6 

100. 

1 . 44  175.8 

899.3  .0  2388.4 

.6  - .4. 

-.0. 

-.03 

5 

760.6 

98. 

.95  177.9 

899.8  .0  6.76.4 

1.4  1.3 

.0 

.03 

5 

770.6 

99. 

„ 77-  1.79.5 

899.8  - -.  .0  -.453. 0 

2.2.  2.0-  _ 

-_0-  03 

5 

target 

IMPACT 

'AT  180.  NM 

AND  TIME  774.6 
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S 
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.<;  s s s 

S ! 

s s 

s s 

SSSS  SSSSS  SSSSS  SSSSS  SSSSS 

s 0 
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0. 

.72  44.9 

899.8  .0  404.1 

2.5  2.2 

.0 

.03 
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